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Part  I 
Liquids  with  Two  Layers  of  Different  Density 

1.   Intro due ti on 

The  main  purpose  of  these  papers  is  to  discuss  what  are  called 
solitary  waves  in  gravitating  incompressible  liquids  of  non-constant 
density.   These  are  steady  two-dimensional  flov7S  with  a  free  surface 
in  a  channel  with  a  horizontal  bottora  vjhich  extends  to  infinity  in 
both  directions.   The  essential  characteristic  of  such  waves  is  that 
they  have  a  sin^jle  crest  vihile  the  vertical  displacement  of  the  free 
surface  tends  to  the  equilibrium  level  at  infinity.   For  liquids  of 
constant  density  the  problem  has  a  lengthy  history  which  will  be 
outlined  a  little  later  in  this  introduction. 

Before  recounting  this  history,  ho^^rever,  it  is  of  interest  to 
examine  the  problem  in  a  general  way  with  reference  to  the  underlying 
mathematical  theory;  first  of  all  for  the  kno^'in  case  of  a  single 
layer  having  constant  density,   A  velocity  potential  ^(x,y)  is  to  be 
found  for  the  space  occupied  by  the  liquid.  As  is  well  known,  this 
function  is  harmonic, 


(1.1)  ^   +  (i)   =  0  , 
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and  it  satisfies  the  condition 


(1.2)  L(x»0)  =  0 

at  the  bottom  of  the  channel.   At  the  free  surface,  characterized  by 
the  equation  y  =  y  (x),  (see  Pi^^.  1.1)  the  boundary  conditions  are 


J  =ys(-) 
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Pic.  1.1 


(1.3) 
(1.1+) 


^y^-ys^-^sA^^'J^s^  =  ° 


2sy3  +  l^(x,yg)+5y(x,y^)  =  2gh  +  c^  . 


The  first  of  these  conditions  is  the  kineraatic  condition,  the  second 
is  Bernoulli's  law,  with  the  pressure  presumed  zero  at  the  free  sur- 
face.  At  infinity  the  conditions  are  that  y  — >  h  and  "^  — >  c,  a 
constant.   Even  though  our  main  purpose  is  to  discuss  solitary  waves, 
it  has  some  point  to  discuss  here  also  the  case  of  vjaves  xvhich  are 
periodic  in  x. 

Because  of  the  unknown  position  of  the  free  surface  and  the  non- 
linear boundary  conditions  there,  no  frontal  attack  on  the  problem 
formulated  above  has  been  successful.  Hox\rever,  procedures  involving 
developments  with  respect  to  a  small  parameter  have  been  employed  to 
obtain  flov;s  in  a  neighborhood  of  the  flows  with  constant  horizontal 
velocity  everyvjhere,  and  a  horizontal  free  surface  y  =  h.   These 
constant  flows  are,  as  one  sees  imined lately,  exact  solutions  of  the 


nonlinear  problem:  "^(x,y)  =  ex,  7_(x)  =  h  satisfy  all  the  conditions 
(1,1)  to  (1,!;),  and  the  conditions  at  infinity.   If  the  depth  h  is 
fixed,  we  have  therefore  a  one -parameter  family  of  exact  solutions 
of  the  problem  with  the  velocity  c  as  paraiaeter. 

Two  quite  different  procedures  have  been  used  to  find  solutions 
of  the  problem  in  the  neighborhood  of  these  flows.   The  first,  and 
more  direct,  of  these  is  a  perturbation  procedure  with  respect  to  the 
amplitude  of  the  disturbance;  this  theory  thus  Is  in  the  category  of 
small  vibration  theory  of  a  continuum  near  an  equilibrium  position. 
One  simply  assumes  for  ^  and  y  developments  of  the  form 

(1.5)  l(x,y)  =  CX+ e42(x»y)  +  e  k^Sy^il)-^   ...   , 

(1.6)  yg(x)  =  h+ Ey^(x)  +  eS^^^^)"*"  •••   » 

and  inserts  them  in  all  of  the  basic  equations,  iN^hich  then  furnish 
a  sequence  of  linear  problems  for  the  successive  determination  of  the 
coefficients  of  the  series.   A  consequence  of  this  procedure  is  that 
the  free  surface  conditions  are  to  be  satisfied  at  y  =  h,  that  is, 
at  the  equilibrium  position  of  the  free  surface  (as  always  in  prob- 
lems of  small  vibrations  about  an  equilibriujn  position).   For  the 
lowest  order  terms  the  free  surface  conditions  are 

(1.7)  4iy(x,h)  -cy^^  --=  0  , 

(1.8)  gy^  +c4-Ly(x,h)  =  0  . 
The  elimination  of  y-,  (x)  gives 

(1.9)  g4i^(x,h) +c^4i,^x^^»^^  ""  °  • 


Of  course,  all  functions  (\.{x,j)   are  hsjnrnonic  functions  In  the  strip 
-co  <x<oo,  0<y<h  subject  to  c^^^^C^jO)  =  0.   Both  the  4-  and  y. 
are  assumed  to  be  bounded  at  infinity.   There  exist  periodic  solutions 
of  this  homogeneous  linear  problem  for  4-1  (x,y),  aside  from  the 
"trivial"  solutions  4i  =  0,  y,  =  0,  namely 

(1.10)  4n (x,y)  =  A  cosh  my  •  sin  m(x+a)  , 

with  A,  a  arbitrary'-  constants,  provided  that 

/T  T  T  N                  tanh  mh  _  c  ' 
U.ll)  — ^^1^; -,  ^-  . 

These  solutions  have  been  shoMi  by  VJeinstein  [22]  to  be  the  only 
solutions  of  the  problem  vjhich  are  bounded  at  infinity.   The  condition 

(1.11)  insures  that  the  free  surface  condition  (1.9)  is  satisfied. 
Once  c  and  h  are  c^-ven,  the  equation  (1,11)  determines  the  wave 
number  m  or  the  wave  length  T  from  T  =  27i/m  .   The  equation  (1.11) 
also  displays  the  important  dimensionless  parameter  c  /gh,  which 
fixes  the  character  of  the  solutions.   The  solutions  are  quite  dif- 
ferent  depending  upon  whether  c  /gh  is  less  than,  equal  to,  or 
greater  than  1;  and  in  fact  the  situation  with  respect  to  the  linear 
problem  is  well  known  to  be  as  follows  (cf.  for  example  Stoker  [2li.] , 
p.  207); 

1.  c  /gh  <  1  .   A  oer iodic  wave  exists.   For  fixed  h  the  wave- 

o 

length  T  Increases  ixrith  c  and  tends  to  co  as  c'^/gh  — >  1. 

2,  cvgh  >  1  .   Only  the  trivial  solution  exists. 

In  the  first  case  it  har.  been  proved  by  Ctruik  [23],  following 
the  method  of  Lcvi-Civita  [10]  for  infinite  depth,  that  rigorous 
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periodic  solutions  of  the  exact  nonlinear  problem  exist  by  showing, 
in  effect,  that  the  series  (1.5)  and  (1.6)  converge  for  e  sufficiently 
small.   Thus  the  waves  furnished  by  the  linear  theorj''  are  indeed  the 
lowest  order  approximation  in  a  convergent  development  with  respect 
to  amplitude.   There  are  bifurcations  of  the  solutions  at  every  value 
of  c  with  c  /gh  <  1.   Indeed,  more  than  one  exact  solution  exists  for 
these  parameter  values.   In  the  second  case  no  solutions  of  the  non- 
linear problem  except  the  trivial  solution  can  be  found  by  this 
method — all  of  the  4-  and  y.  would  vanish. 

It  is  worthwhile  to  state  what  happens  if  the  problem  for  4t 
and  y,  is  made  nonhomogeneous,  say  by  inserting  a  small  obstacle  in 
the  bed  of  the  stream,  or  by  applying  a  small  nonzero  pressure  over 
a  portion  of  the  free  surface  to  create  a  disturbance.   In  case  1) 
no  unique  solution  can  exist,  clearly,  since  the  solution  (1.10)  of 
the  homogeneous  problem  can  always  be  addedj  in  this  case  a  unique 
solution  is  picked  out  in  practice  by  assuming  a  radiation  condition 
at  infinity,  i.e.,  that  the  disturbance  created  dies  out  far  upstream. 
In  case  2),  no  radiation  condition  need  be  applied,  since  the  homo- 
geneous problem  has  only  the  trivial  solution.   However,  in  the 
special  case  c  /gh  =  1,  the  steady  flow  solution  for  4-i  becomes  un- 
bounded at  infinity  in  contradiction  with  the  linearization.   Indeed, 
as  we  shall  have  ample  occasion  to  see,  flows  at  and  near  the 
critical  speed  given  by  c  =  s,/gH  ,  which  is  the  upper  bound  for  the 
speed  of  sinusoidal  progressing  waves  In  still  xv'ater,  are  special  in 
character  and  require  a  special  treatment.   However,  L.  Nlrenberg  has 
proved,  in  an  unpublished  paper,  that  rigorous  solutions  of  the  non- 
linear problem  exist  both  if  c  < s/gh,  and  if  c  > Vgh  when  obstacles 
of  small  amplitude  (i.e.,  of  first  order  in  e)  are  put  in  the  stream. 
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In  addition  to  the  steady  flows  it  is  also  interesting  to  con- 
sider the  linear  theory  of  the  unsteady  flows  which  result  v/hen  small 
obstacles  are  put  in  the  streari  vjith  uniform  flow  at  t  =  0  and  it  is 
required  to  describe  how  the  disturbance  builds  up  as  time  goes  on. 
If  0  '/gh  r   1>  o^G  finds  as  t  -?>  oo  the  expected  result  that  the 
solutions  tend  to  the  steady  flows  which  were  in  question  in  the 
preceding  paragraphs.   However,  if  c  /gh  =  1,  i.e.,  if  the  undisturbed 
flovj  has  the  critical  speed,  It  is  found  (cf.  [2^4.])  that  the  linear 
theory  yields  a  disturbance  i-^hlch  becomes  large  everywhere  without 
limit  as  t  — >  00  ,   One  would  be  tempted  to  say  that  the  flow  with 
critical  speed  is  very  unstable--the  procedure  formulated  by  the 
linear  theory  is,  in  fact,  the  standard  small  variation  method  for 
discussing  the  stability  of  a  motion.   The  small  oscillation  crite- 
rion fails,  however,  in  this  case,  it  would  seem,  from  both  the 
mathematical  and  exoerlmental  points  of  view,  since  the  vjave  called 
the  solitary  vjave  exists  in  the  vicinity  of  the  critical  speed.  Prom 
the  point  of  view  of  observations  and  of  laboratory  experiments,  this 
wave  seems  to  be  extremely  stable,  and  it  also  seems  rather  easy  to 
create  it. 

Recently,  Priedrichs  and  Ilyers  [16]  hPve  proved  rigorously  the 
mathematical  existence  of  a  solitary  vjave  solution  of  the  full  non- 
linear problem.   The  method  they  use  has  been  referred  to  above  as 
the  second  of  tvjo  basic  procedures.   It  is  not  a  method  of  perturba- 
tion with  respect  to  amplitude.   Instead,  it  is  an  iteration  pro- 
cedure based  on  the  small  parameter  o"  =  c  /gh  -1,  and  it  therefore 
furnishes  solutions  valid  near  the  uniform  flow  at  critical  speed. 
This  procedure,  hovjever,  requires  that  the  dependent  and  Independent 


variables  be  stretched  In  viays  i.'hlch  depend  oner,  and  also  in  a 
manner  which  distinguishes  the  horizontal  frora  the  vertical  direction. 
To  lowest  order,  however,  it  furniahej.'  the  sarae  result  as  a  formal 
perturbation  procedure  obtained  by  setting  up  power  series  in  CT' 
(which  are  set  up  after  the  stretching  process,  of  course);  this  was 
carried  out  earlier  by  Keller [13].   In  fact,  to  lowest  order,  the 
theory  is  identical  with  the  approximate  theory  called  the  nonlinear 
long  wave  or  shallow  water  theory.   The  result  in  the  lowest  order 
of  approximation  is  essentially  the  s.^ine  as  the  well-knov;n  much 
earlier  results  of  Boussinesq  [3],  Rayleigh  ['4],  and  others. 

It  is  highly  unlikely  that  perturbation  series  in  powers  of  the 
stretching  parameter  O^"  would  converge;  such  series  seem  certain  to 
have  at  best  asymptotic  significance  since  the  stretching  process  has 
the  effect  of  making  X  a  coefficient  of  some  of  the  highest  deriva- 
tives occurring  in  the  problem,  so  that  boundary  layer  effects  are 
to  be  expected  in  the  limit,  and  the  perturbation  procedure  seems 
not  likely  to  be  of  the  regular  type.   In  this  paper,  formal  series 
in  the  parameter  cr  are  used  to  obtain  the  lowest  order  terms. 

In  this  manner  one  can  find  solitary  waves  of  arbitrarily  small 
amplitude;  but  these  vmves  cannot  be  obtained  by  perturbation  series 
in  the  amplitude  in  the  way  described  above,  as  was  already  remarked: 
the  coefficients  of  such  series,  on  the  basis  of  the  uniqueness 
theorem  of  VJeinstein,  would  all  vanish,   Thu.^.  the  solitary  wave  is 
not  analytic  in  the  amplitude. 

Although  only  passing  remarks  are  made  about  so  called  cnoidal 
wave  solutions  in  this  paper,  it  is  perhar)s  vjorthwhile  to  mention 
them  briefly  here.   In  addition  to  the  solitarj'^  vjave,  there  exist 
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also  periodic  waves  in  the  neighborhood  of  the  critical  speed.   These 
waves  are  obtained  by  the  same  means,  described  above,  which  lead  to 
the  solitary  i\fave ,  and  Littman  [lO]  has  proved  the  existence  of  such 
waves  by  using  a  method  similar  to  that  of  Friedrichs  and  Hyers.   For 
small  departures  from  the  critical  speed  these  waves  have  large  wave 
lengths  and  in  the  limit  as  the  wave  length  goes  to  infinity  they 
tend  to  the  solitary  wave. 

The  object  of  this  paper  is  to  treat  problems  of  the  same  gener- 
al types  as  those  first  discussed,  but  for  liquids  of  non-constant 
density.   In  doing  so  attention  is  centered  mostly  on  solutions  of 
the  type  of  the  solitary  v.'ave  or  cnoidal  wave.   I^o  convergence  proofs 
are  attempted;  lowest  order  terras  only  in  formal  perturbation  schemes 
are  obtained.   However,  it  seems  highly  probable  that  these  solutions 
are  indeed  the  correct  lowest  order  terms  in  the  development  of  exact 
solutions  of  the  nonlinear  problem  with  respect  to  the  appropriate 
parameters . 

In  Part  I,  the  case  of  a  fluid  vjith  two  layers  each  of  constant 
density,  vdth  a  free  surface  on  top,  is  treated.   Both  the  linear 
theory  of  waves  of  small  amplitude,  and  a  nonlinear  theory  of  the 
same  kind  as  the  long  wave  or  shallow  water  theory  are  employed.   The 
linear  theory  furnishes  the  important  result  in  this  case  that  there 
are  two  critical  speeds.   These  speeds  are  defined  as  the  limiting 
speeds  of  waves  whose  vjave  length  tends  to  Infinity,   For  layers  of 
equal  thickness  h/2,  for  example,  but  different  densities  6,  /\   such 
that  if  5  is  the  density  of  the  lox>.'est  layer  p  =  A,/5  <  1,  these 
speeds  are  given  by 
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Thus  if  the  densities  sre  nearly  the  saviie  the  larger  critical  speed 
is  near >/gh  as  one  might  expect,  and  the  smaller  is  near  zero.  From 
the  discussion  above  one  might  also  exoect  that  two  solitary  vjaves 
would  exist  with  speeds  near  the  critical  speeds.   That  turns  out  to 
be  the  case.   The  tv/o  critical  speeds  distinguish  tx-jo  types  of  soli- 
tary waves.   One  type  has  characteristics  similar  to  the  ordinary 
solitary  wave.   The  stream  lines  for  this  type  are  lines  of  elevation 
above  the  undisturbed  stream  lines  and  they  have  a  maximum  amplitude 
at  the  free  surface.   The  other  type  tends  to  be  hidden  from  observa- 
tion at  the  free  surface.   The  waves  are  internal  with  a  maximum 
amplitude  at  the  interface  much  larger  than  the  amplitude  at  the  free 
surface.   Depending  on  a  value  of  a  function  of  the  ratio  of  the 
depths  and  the  ratio  of  the  densities;  of  the  layers,  the  internal 
vjave  at  the  interface  may  be  either  a  vjave  of  depression  or  elevation, 
Cnoidal  waves  near  the  two  critical  speeds  can  also  be  determined. 
In  Part  II,  where  fluids  vjith  a  vsrtical  densitjr  distribution 
which  decreases  exponentially  going  upward  are  treated,  it  is  found 
that  the  linear  theory  for  waves  of  small  ajnplltude  furnishes  an 
infinite  spectrum  of  critical  speeds.   These  speeds  form  a  discrete 
set  with  0  as  the  only  limit  point.   One  is  led  to  conjecture  that 
there  would  be  infinitely  many  solitary  (and  also  cnoidal)  waves 
associated  with  the  critical  speeds  and  that  the  extrema  of  each  of 
these  waves  would  be  associated  with  different  levels  in  the  fluid 
which  would  increase  in  nuraber  and  begin  nearer  the  bottom  as  the 
corresponding  critical  speed  Is  taken  siaaller.   These  conjectures 
all  prove  to  be  true. 
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We  turn  now  to  a  discussion  of  the  earlier  history  of  this  sub- 
ject.  The  main  object  of  the  account  which  follows  is  to  briefly 
indicate  the  basic  methods  vrhich  hcve  been  used  in  the  developraent 
of  the  theory  of  solitary  and  cnoidal  vjaves,  and  to  point  to  the 
fundaraental  theoretical  papers  ^^Ihich  to  our  knowledse  are  concerned 
primarily  with  solitary  or  cnoidal  ^^javes.   There  are  papers  vjhich 
are  concerned  with  the  application  of  the  solitary  wave  theory  to 
ocean  waves  and  atmospheric  vjaves,  and  there  are  others  x-jhich  report 
results  of  laboratory  experiments;  some  of  these  are  listed  tov;ard 
the  end  of  the  bibliography. 

The  name  solitary  vjave  is  due  to  Scott  Russell  [1]  xi/ho  presented 
the  first  scientific  report  (iSlili)  about  such  a  vjave.   Scott  Russell 
vjas  an  engineer  who  studied  the  solitary  iNfave  in  connection  with  an 
investigation  designed  to  ascertain  the  feasibility  of  operating 
steam  driven  craft  in  the  Edinburgh  and  Glasgox-j  canal.   His  observa- 
tions led  Russell  to  adopt  the  empirical  formula 

c  =yg(h+aT 

for  the  velocity  c  of  the  solitary  wave  in  terms  of  g,  the  accelera- 
tion due  to  gravity;  h  the  undisturbed  depth  of  the  channel;  and  a, 
the  am.plitude  of  the  wave  T\'hich  Is  not  necessarily  small.   This 
formula  was  confirmed  by  the  more  lorecise  experiments  conducted  by 
Bazln  [2]  in  1865.   The  first  mathematical  theories  which  gave 
approximate  descriptions  of  the  profile  and  speed  of  the  solitary 
wave  vjere  presented  independently  by  Boussinesq  [3]  (I67I)  and 
Rayleigh  [l\.]    (I876).  Rayleigh's  analysis  proceeds  from  the  expansion 
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of  a  tirae  independent  coraplex  velocity  potential  in  powers  of  the 
vertical  coordinate.   Boussinesq  used  the  first  terras  In  a  similar 
expansion  of  a  tirae  dependent  velocity  potential  and  then  supple- 
mented his  analysis  with  various  physical  assuraptions  in  addition 
to  those  iraplied  by  the  basic  hydrodynai'nlcal  theory, 

McCowan  [5]  (1891)  introduced  a  coraplex  velocity  potential  vjhich 
gave  a  second  order  approximation  to  the  speed  of  the  solitary  wave. 
Kortei^'eg  and  DeVries  [6]  (1895)  used  a  raodif ication  of  Rayleigh's 
method  to  shovj  that  a  general  class  of  long  waves  of  finite  amplitude 
and  permanent  type  raay  exist.   They  found  waves  described  by  the 
Jacobian  elliptic  function  cn(^,k).   The  existence  of  these  periodic 
waves  which  travel  without  change  of  form  was  first  indicated  by 
Boussinesq  [?]  (I89I).   Korteweg  and  DeVries  called  these  waves 
cnoidal  waves  and  showed  that  the  limit  of  a  cnoidal  wave,  whose 
wave  length  approaches  infinity,  is  a  solitary  wave.   Gwyther  [8] 
(1900)  and  Levi-Civita  [9]  (190?)  presented  refinements  and  system- 
izations  of  Rayleigh's  method. 

In  1925  Levi-Civita  [10]  proved  the  mathematical  existence  of 
two-dimensional  periodic  waves  of  finite  amplitude  and  permanent 
type  in  water  of  infinite  depth.   Levi-Civita  used  the  velocity 
potential  function  and  the  stream  function  as  independent  variables. 
With  these  variables  the  problera  was  reduced  to  finding  a  function 
analytic  in  a  known  doraain  and  subject  to  a  certain  boundary  condi- 
tion nonlinear  over  part  of  the  boundary.  Weinsteln  [11]  (1926) 
used  Levi-Civita' s  raethod  to  improve  the  known  approximations  to  the 
speed  of  the  solitary''  wave.   A  numerical  mistake  in  Weinstein's 
second  order  forraula  has  been  corrected  by  Long  [20].   Davies  [12] 
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(1951)  investigated  waves  of  finite  craplitude  by  modifying  the  non- 
linear condition  vjhich  appears  in  the  Levi-Civita  presentation;  and 
then  Packham  [13]  (1952)j  using  Davies'  technique,  studied  the 
solitary  wave  and  found  a  forniula  for  the  speed  similar  to  the  one 
found  by  licCowan. 

In  an  aiopendix  to  a  psoer  by  Stoker  [llj.]  (19^3)  Priedrichs  gave 
a  systematic  perturbation  raethod  designed  to  produce  the  conventional 
non-linear  shallow  vjater  theor^'-  in  the  first  order  approximation,  and 
to  exhibit  a  parameter  which  determines  the  accuracy  of  this  and 
higher  order  approximations.   Friedrichs'  method  is  characterized 
by  the  introduction  of  a  dimensionless  parameter  cr~  which  he  used  to 
stretch  some  of  the  variables;  and  the  assumption  that  all  quantities 
can  be  expanded  in  powers  of  a-.   Keller  [l5]  (19lj.8)  used  Friedrichs' 
method  to  show  that  in  the  case  of  steady  motion,  the  second  order 
aporoximation  gives  cnoidal  and  solitary  waves.   Keller  developed  a 
systematical  procedure  for  the  theoretical  study  of  these  waves  and 
found  results  in  essential  agreement  with  those  found  by  Rayleigh, 
Boussinesq,  and  Korteweg  and  DeVries.   Later  Priedrichs  and  Hj'-ers 
[16]  (1951|)  used  Friedrichs'  method  to  prove  the  matherr^atical  existence 
of  the  solitary  wave.   In  [1&]  Lavrentiev  [1?]  (19'46),  who  used  the 
Levi-Civita  approach,  is  credited  as  being  the  first  to  give  a  proof 
of  the  existence  of  the  solitary  wave.   Littman  [18]  (1957)  used  a 
method  similar  to  the  one  in  [16]  to  prove  the  existence  of  cnoidal 
waves  near  critical  speed. 

In  each  of  the  papers  mentioned  above  it  was  assumed  that  the 
density  of  the  medium  is  constant.   In  recent  years,  there  has  been 
a  growing  interest  in  solitary  v/aves  in  a  stratified  medium.   The 
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interest  stems  frora  the  surmise  that  solitary  waves  may  play  a  con- 
siderable role  in  certain  meteorological  "ohenomena  and  phenomena  in 
other  fields  vjhere  density  changes  are  important.   Keulegan  [19] 
(1953)  used  Boussineeq's  method  to  investigate  the  characteristics 
of  internal  solitary  waves  at  the  interface  of  a  txra-layer  medium 
bounded  by  rigid  planes  at  both  the  top  and  bottom  surfaces.   Long 
[20]  (1956)  used  Rayleigh's  method  to  find  a  fourth  order  approxima- 
tion to  the  speed  of  the  ordinary  solitarj'-  wave,  and  in  the  sarae 
paper  he  verifies  some  of  Keulegan' s  results  and  gives  an  second 
order  approximation  to  the  internal  solitary  wave  of  a  two-fluid 
system,  again  bounded  by  rigid  planes  above  and  belovj.   Abdullah  [21] 
(1955)  used  Keller's  procedure  to  discuss  the  solitary  wave  at  the 
interface  of  a  two-layer  medium  for  which  it  was  assumed  that  the 
hydrostatic  law  holds  in  the  upper  layer. 

In  the  present  part  a  variant  of  Priedrichs'  method  is  used  to 
obtain  solitary  waves  in  a  two-layer  nsdium  with  a  free  surface  at 
the  top  of  the  upper  layer.   In  Part  II  v.'e  will  present  results 
which  have  been  found  for  solitary  waves  in  a  medium  in  which  there 
is  a  continuous  variation  of  densit^r.   In  addition  to  results  already 
mentioned,  we  present  in  Part  I  fori.iulas  for  the  velocities  of  the 
two  kinds  of  waves  in  terms  of  the  maximum  amplitudes,  and  expressions 
for  the  streajmlines  throughout  the  medium  as  vjell  as  expressions  for 
the  velocity  components  and  pressure  of  the  flow. 


i.I 


Ill 


Por:'"ulation  of  the  Problem 


In  order  to  study  tvjo-dlraenslono.!  x-javes  of  permanent  type  in  a 
two-layer  mediura  of  iriuniscible  liquids  let  us  suppose  that  a  cross 
section  of  the  medium  in  the  equilibrium  state  fills  an  infinite 
horizontal  strip.   The  lower  layer  whose  density  is  constant  and 
equal  to  S  is  supposed  to  be  supported  by  a  rigid  bottom;  and  the 
upper  layer,  x-jhose  density  is  constant  and  equal  to  ^  such  that 
/\^  <  5,  is  supposed  to  have  a  free  surface,  i.e.,  a  surface  on  which 
the  pressure  p  is  zero  and  there  are  no  geometric  constraints.   Let 
us  supoose  that  a  disturbance  of  some  kind  in  the  medium  initially 
at  rest  has  created  a  wave  of  -oermansnt  type  which  moves  to  the  left 
with  velocity  c;  and  that  the  character  of  this  wave  is  observed 
from  a  coordinate  system  which  moves  with  the  wave.   The  x-axis  of 
this  system  is  taken  to  coincide  with  the  bottom,  and  the  vertical 
y-axis  is  chosen  so  as  to  pass  through  a  crest  of  the  wave.   VJith 
respect  to  these  axes  the  wave  is  stationary  and  the  velocity  of  the 
medium  at  infinity  is  c.   As  shown  in  the  Pig.  2.1,  let  us  take  the 
depths,  or  the  mean  depths,  of  the  lov;er  and  upper  layers  at  infinity 
to  be  h  and  H  respectively. 


/  / 


/   /    /  /    ■'■  .'  /  /    ^  / 


/  / 


/  /   /  /    /    /  /  / 


X 


Fis.  2-.1 
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The  steady  state  equations  vrhich  have  to  be  satisfied  for  the 
lovjer  layer  are 

0 ( UU  +VU  )  =  -p   , 

5(uv^+vVy)  =  -og  -Py.  , 


/  w    /v 


u  +v   =  0  , 
X  y     ' 


/■^'V'      /^.^/*-' 


US  +VS   =  0  , 
X   y 


-ij 


^  =  V^-^y  > 
where  u,v  are  the  velocity  components;  p  is  the  pressure;  g  is  the 
acceleration  due  to  gravity;  s(x,y)  =  0  is  the  equation  of  a  boundary 
curve;  and  ^  is  the  vorticity.   The  doraain  in  vjhich  the  solution  is 
to  be  found  is  partly  unknown  at  the  outset.   If  the  density  6  is 
replaced  by  /^  v;e  have  the  equations  for  the  upper  layer.   Hereafter, 
quantities  which  pertain  to  the  upoer  layer,  and  v/hich  need  to  be 
distinguished  from  similar  quantities  in  the  lov;er  layer,  will  be 
denoted  by  capital  letters. 

If  the  stream  function  ^'   is  \{/(x,y)  so  that 

u*  =  i!/y{x,y)  V  =  -\!/^(x,y) 

the  streraa  lines  are  given  implicitly  by 

i!/(x,y)  =  Y  =  const. 
The  function  \!/(x,y)  can  be  chosen  so  that  for  the  bottom 

Y  =  i!/(x,0)  =  0 


(-) 


The  streara  function,  rather  than  the  velocity  potential,  is 
Introduced  because  it  is  convenient  here,  and  more  or  less  in- 
dispensable for  the  case  of  continuous  variation  in  density  to 
be  treated  in  Part  II. 
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and  then  the  value  of  y  which  determinss  the  interface  can  be  taken 
to  be 

Y  =  ch  , 

while  for  the  free  surface  we  rnay  write 

Y  =  c(h+H)  =  ch(l+r) 

v/here  r  =  H/h.   It  is  assumed  that  to  each  value  of  y  there  corre- 
sponds one  and  only  one  stream  line  idiich  is  rriven  by  solvin,'^ 
Y  =  \!/(x,y)  for  y  namely, 

y  =  r(:c,Y)  . 

The  function  f (x,y)  will  be  called  the  stream  line  function.   This 
assumption  is  reasonable  since  flows  near  uniform  flovjs  are  the  only 
ones  considered. 

Let  us  introduce  x  and  y  ^s  independent  variables  and  use  the 
notation 

u(x,y)  =  u[x,f(x,Y)]  =  u(x,y) 

and  a  similar  bar  notation  for  the  other  transformed  circumflex 
quantities.   Then,  for  example 

u^(x,y)  =  u^+u^  11  =  u^-v5^ 

and  using  these  relations  the  momentum  equations  for  the  lower  layer 
become 

5uu^  =  -Px+vPy  , 
o^v^  =  -5g  -  up;^.  , 
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the  continuity  equation  becomes 

u  -vu  +UV   =  0  , 
X   Y   Y 

and  the  vorticity  is  given  by 

The  kinematic  condition  at  a  boundary  curve  becomes 

u(s^-vs^)+vuSy  =  0 
or 


us   =  0 

X 

but  this  of  course  is  automatically  satisfied  since  for  a  boundary 
curve  we  have  's(x,y)  -   0  where  y  is  a  constant  corresponding  to  the 
boundary  stream  line.   The  domain  for  xjhich  the  differential  equations 
must  be  satisfied  is  now  the  strip  -co  <x<oo;  0<y<  ch.   By 
working  with  the  stream  function  in  this  wa^r  the  advantage  of  a  fixed 
domain  for  the  independent  variables  is  achieved. 

Instead  of  using  the  dependent  variables  u,  v  and  p  it  is  more 
convenient  for  our  purposes  to  use  u,  f  and  p  as  the  dependent 
variables.   For  these  variables  we  have 

7  =  -L  =  1  r     =   "  X  _  V 

Y   ^    -  X    ij/     - 

'    ^y   u  T    u 

and  it  can  easily  be  verified  that  the  continuity  equation 


u  -vu  +UV   =  0 
X   Y   Y 
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is  automatically  satisfied,   ^'"e  do  not  expect  u  to  vanish  since  the 
flow  is  near  the  uniform  flovj  uith  velocity  c  ^  0.   The  momentum 
equations  nov;  take  the  form 

5u  =  -f  p  +f  p 
X     Y  X   x^  Y 

5(uf   +u  f  )  =  -5gf  -  p  , 

*      XX       XX  O   y     fy 

and  with  these  i-ie   have 

Uf y     =     1     , 

as  a  replacement  for  the  continuity  equation.   The  vorticity  is 

given  by 

—  —2 

f     ,  ^   1+f 

Y  ^y 

Finally,  if  we  introduce  the  dimenslonless  variables 


«  =  i 

'>  =  h 

---c 

"  =  l 

^  =  1 

v  =  ! 

p  =  £-^ 

'  -  Ae^ 

^  =  1 

^  =  1 

P  =  # 

^  = 

_  h 
c 

K 

we  find  that  the  differential  equations  become 


(2.1)  uf^^  +  u^f^  =  -  Xf^  -  p^^  , 

uf^  =  1  . 
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They  hold  in  the  strip  -co  <x<oo;0</;f<l.   In  the  upper  layer 
the  equations  are,  in  the  same  v/ay; 


"5  =  -  \h  *  "^n 


(2.2)  UP 


UP^  =  1  . 


and   they  hold   in  the   strip   -00    <x<co;l<1<  1+r.      The  boundary 
condition  at   the  bottom   is 

f(x,0)    =  0    . 

The  boundary  conditions  at  the  interface  are 

f(x,l)  =  P(x,l)  ;    p(x,l)  =  pp(x,l)  . 

These  conditions  state  that  the  interface  is  a  stream  line  and  that 
the  pressure  is  continuous  along  it.   Since  the  pressure  at  the  free 
surface  is  zero  the  boundary  condition  there  is 

P(x,l+r)  =  0  . 

For  irrotational  motion,  vjhich  we  assume,  the  vortlcity  must  be  zero 
in  each  layer  and  this  means  that  both  f  and  P  must  satisfy 

^    f.       n   3    1+f? 

or  ^ 

(2.3)  (l+f?)f.^  +  fff^^  -2f^f^f.„  =  0  . 

The  equation  ^^7hich  expresses  Bernoulli's  law  is  found  by  elim- 
inating the  horizontal  velocity  component  from  equations  (2.1),  (2.2) 
and  Intesratingj  for  example: 


i^  >   I  X  >    oo-   rJ^'i'^'"-^ 


<    X    W   *^  Ck 
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•3  ,.1_ 


af  (f7^  =-  V^'V^ 


I  ( — ^)  +  p  +  xf  =  an)  . 

For  irrotational  flow  C  (^' )  is  a  constant. 


3«   The  Linearized  Equations.   Critical  Speeds 

For  the  case  of  a  steady  parallel  flow  of  velocity  c  in  each 
layer  we  have 

u  =  1  U  =  1 

f  =  /?  F  =  ^ 

p  =  \{l+pr'{)  P  =  X(l+r-^J  . 

We  consider  the  vjave  motion  to  be  a  sraall  disturbance  superposed  on 
this  state  of  constant  flow.   We  therefore  i^rrite 

u  =  1+u'"'  U  =  1+U'"" 

f  =  n+f'"'  F  =  '2+F" 

p  =  X(l+pr-'Z)+p"''       P  =  7^(l+r-1)+P''''  . 

Then  if  we  substitute  these  quantities  in  (2.1),  (2.2),  (2,3)  and 
assume  the  disturbance  so  sraall  that  second  order  products  involving 
the  starred  quantities  can  be  neglected  vie   find  the  linear  equations 


»v;: 
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u"+f,;  =  0 


subject  to  the  boundary  conditions 


u; 


-P-.XP^ 


f""'(?,0)  =  0 

f'''(^,l)  =  p'""(5,l)  p''(^,l)  =  pp'''"(c,l) 

P'"(^,l+r)    =   0    . 
For   irrotational  motion  f  and  P  raust    satisfy 


(3.1) 


f       +  f 


^k  -^  ^^^ 


=   0 


=   0 


•00     <   ^    <  CO 

0  <    *l  <  1 

•00     <    5    <  00 

1  <  ^  <  1+r 


and   the  bouiidary  conditions 


.--^ 


f"(^,l)   =  P"(^,l) 


f'""(5,0)    =  0 

Xf|(C,l)-fJ^(^,l)  =  p[af|(^,i)-p*,^(c,i)] 
F;f^(C,l+r)   =   XP|(^,l+r)    . 


The  following  functions  are  solutions  of  (3.1)  which  satisfy  the 
boundary  conditions  at  ''(   =0  and  ^=1: 


f'"'  =  A  sin  (v^+B)  sinh  v 


(3.2) 


P  = 


A_ 
pv 


sin  (vC+B) 


V  sinh  v^-v(l-p)  sinh  v  cosh  v(l-'t) 
+  A(l-p)  sinh  V  sinh  v(l-''/) 


with  V,  A,  and  B  any  constants. 


-     V 


»v  rlShi 


There  are  solutions  which  correspond  to  the  special  value  v  =  0 
of  V  and  which  satisfy  the  boundary  conditions  at    fl  =  0   and  ^  =  1. 
They  are : 

f  =  {at,+b)f-i 

F  =  h^±l    [(l-A+XpR+(l-p)(X-l)]  . 

If  a  is  not  zero,  these  solutions  .rrive  stream  lines  which  depart  as 
far  as  v.-e  please  from  f  =^L;   F  =^  and  therefore  they  cannot  be 
admitted  in  the  linear  theorj'-. 

If  the  functions  in  (3.2)  are  to  satisfy  the  boundary  condition 
at  'l  =  1+r  for  v  ^  0  one  must  have  (1+p  tanh  v  tanh  vr)v 
(tanh  v+  tanh  vr)\v+X    (1-p)  tanh  v  tanh  vr  =  0  which  can  be  written 

(3.3)     (1+p  tanh  V  tanh  vr)[^  -  a^(v)][-^  -  a^{v)]      =  0 


where 


/(tanh  V  -  tanh  vr ) 

(tanhv+tanhvr )+  !_.•  .  ,  ■  ,  n  /•,  \*.  ,.4.  i„  i 
„  ,  \  „  ,  \  _  -  \/ +Liptanhvcanhvr  [l-(l-p)tanhvtanhvrj 
a^(v;,a2(v; 2(l+p  tahh  v  tanh  vr) 


Equation  (3.3)  determines  the  wave  nujnber  in  terms  of  the  ratio  of 
the  depths  of  the  layers  r  =  H/h,  the  ratio  of  the  densities  p  =  /\/5 
and  the  speed,  since  \  =   gh/c  .   If  T  is  the  xvave  length  of  the 
sinusoidal  waves,  v  =  27i/T;  it  is  seen  that  (3.3)  gives  two  speeds 
associated  with  a  prescribed  wsve  length.   The  li'niting  speeds  cor- 
responding to  waves  whose  wave  lenr^ths  approach  infinitj'- — this  was 
one  of  the  waj'-s  used  to  characterize  critical  speeds  in  the  sin:;,'le 


1/  Cj      '     t 


.     i'  '•:-}'. '-S-S.) 


=    ""^ 


■:;j1   \,'   i'n  !!■.■? 


.-'ns  j) 


layer  case--are  given  in  this  case  by 


(3.1;) 


_l_ 


1 


gh 


"2 

gK 


a^(v) 


p     aAv) 


_   l+r+ V(l-r  )''-«-i;pr 


2 


v->0 
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as  one  readily  sees.  For  still  another  reason,  which  will  appear 
below,  these  speeds  c,  and  Cp  are  called  critical  speeds.   It  can  be 
seen  that  the  critical  speeds  satisfy  the  equation 


(3.5) 


X'^(l-p)r-(l+r)X+l  =  0  . 


Conversely,  if  A  corresponds  to  one  of  the  critical  speeds  it 
is  found  that  the  left-hand  side  of  (3.3)  will  contain  the  factor  v^. 
From  this  we  can  see  that  the  linear  theory  is  inapplicable  to  non- 
homogeneous  probleras  in  which  the  speed  is  critical,  just  as  is 
known  for  the  case  of  one  layer.   To  see  this,  consider  a  pressure 
disturbance  which  moves  with  velocity  c  over  the  free  surface  of  the 
two-layer  medium.   The  solution  for  the  stream  lines  in  the  upper 
layer,  for  example,  is 


--"•  -i^v 


P'""  = 


TIQ 


vsinhv'/+v(p-l  )slnhvcoshv(''(-l) 
+  X(p-1  )slnhvsinhv(^( -1) 


dv 


X^coshvcoshvr  (1+ptanhvtanhvr )  [-^  -  a,  (v)]  [-j-  -  ap(v)] 


where  P""  is  the  Fourier  transform  of  the  pressure  at  the  free  surface; 
and  L  is  a  path  along  the  real  axis  except  for  semi-circles,  which 
avoid  the  zeros  of  the  denominator,  and  are  situated  so  that  the 


'j:;- 


-  X     _     I 


iiriv 


<• 


;'.iiij      QSTj     C'j.'  »   i;J\s:.L 


ij.rlw  sonsii 
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disturbance  dies  out  upstream,  l^ov^   if  A  corresponds  to  one  of  the 
critical  speeds  the  denorainator  in  the  integrand,  which  contains 
the  left-hand  side  of  (3.3),  will  have  a  fourth  order  zero  at  v  =  0, 
while  the  niunerator  has  a  second  order  zero  there.   The  effect  of 
this  is  that  when  the  path  is  deformed  to  the  other  side  of  the  real 
axis  in  order  to  discuss  the  disturbance  far  downstream,  the  contri- 
bution from  the  second  order  pole  at  the  origin  will  give  stream 
lines  which  depart  arbitrarily  far  from  the  lines  of  constant  flow. 
This  kind  of  behavior  is  inadmissible  in  the  linear  theory--the 
linear  theory  fails  at  a  critical  speed,  and  we  are  forced  to  turn 
to  a  nonlinear  formulation.   This  is  all  just  as  it  was  for  one  layer, 
We  are  thus  led  to  eiqpect  that  the  two  critical  speeds  are  speeds 
near  which  solitary  waves  will  occur, 

I}.,  Waves  Near  the  Two  Critical  Speeds 

The  non-linear  equations  (2.1)  and  (2,2)  can  be  treated  in 
several  different  ways.   One  way  is  to  eliminate  the  velocities  and 
pressures  so  as  to  obtain  the  equations  which  f (5,'''/)  and  F(^,'0  must 
satisfy  and  then  deal  directly  with  these  equations  subject  to  the 
boundary  conditions  expressed  in  terms  of  f  and  F  only.   For  irrota- 
tional  motion  these  functions  must  satisfy 


(ll.l)  p         - 

(1+f|)P,^;^+F^^P^^-2F^F,^P^;,^  =  0  ,  1  <  '^^<  1+r  , 
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These  equations  could  be  analyzed  by  Raylelgh's  method  vjhich  involves 
expanding  f  and  P  in  powers  of  the  deoth  variable  V«   Another  way  to 
proceed,  along  the  lines  used  by  PrJedrichs  [li;.],  is  to  stretch  the 
horizontal  variable  in  ('-i.l)  by  introducing  the  new  variable 
ye"  E,  =  CT,    and  then  assurae  that 

CO         ,  _ro         , 


ra  1^^  "•  fe=0 

The  authors  used  both  of  these  methods  in  preliminary  work.  Each 
method,  although  effective,  has  certain  technical  disadvantages. 
The  method  vjhich  follows  seems  to  be  best. 

In  order  to  investigate  a  wave  xjhich  moves  with  a  speed  such 
that  X  =  gh/c   is  near  some  value  /,  vje  insert  /  in  the  equations 
(2,1)  and  (2,2)  so  that  the  difference  X-/  appears.   In  other  words 
consider  the  equations 


Ur  =  -  f/,  P^  +  f ,!rP, 


rS 


^'1 


uf^^  +  u^f^  =  (,(^-X)f^  .  /f^  .  p^         0  <  Aj.  <  1 
uf.^  =  1 


and  the  corresponding  set  for  the  upper  layer.   If  we  vrrite 

/  -  A  =  e 
and  introduce  the  new  variable 

the  equations  become 
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(1^.2)  e(uf^_^^+u^f^J  =  ef,^  " /f>,  "  P.,  0  <  ^<  1 


and 


uf,^  =  1 


Uo-  =  -  P,,  P^H-  F^P^ 


ik»3)  e(UF^_^+U.P^)  =  eP^  -  /P^  -  P^  0  <  >-^<  1+r. 


UF^^  =  1 


Since  the  vorticity  is  assuraed  to  be  zero  we  also  have 

Let  us  assume  now  that  all  of  the  quantities  in  the  new 
variables  can  be  expanded  in  integral  powers  of  e,  that  is,  we  write 

CO    , 

f  (c>T<L,  e)  =  ;>   e  A.(^'z)     u  =  ...     p=... 
k=0  ^^ 


k=0 


,.?0.   V 

F(=''',^t,s)  =  ^ e'^F,  (cr,^)     U  =  ...     P  =  ... 


where  the  k   order  approximations  are  subject  to  the  boundary  condi- 
tions 

^1,(^,0)  =  0  , 
(1^.6)       fi,C^,l)  =  F^i(T,l)    ,        Pk(^,l)  =  P\ir-,1)    , 

Pj^(j-,l+r)  =  0  . 

As  a  consequence  of  the  assumption  (If-. 5)  the  vertical  velocity  in  the 
lower  layer,  for  example,  must  have  the  exoansion 


{': 
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V  =  I  =  _-2s  =  uf^  =  y?  ui:^- 

The  idea  is  to  substitute  these  expansions  in  ([l,2),  {l\.,3) ,    {k»k-) 
and  the  boundary  conditions,  and  equate  coefficients  of  likepoi/ers  cf  e. 
This  will  lead  to  a  sequence  of  problems  which  are  to  be  solved.   The 
expectation  is  that  the  analysis  will  show  automatically  that  /  must 
correspond  to  one  of  the  critical  speeds,  and  thus  be  equal  to  either 
X-,  or  7\.„. 

VJe  observe  also  that  the  parameter  e  multiplies  highest  order 
derivatives  in  various  equations,  thus  indicating  that  the  perturba- 
tion is  not  of  the  regular  type.   Rather,  an  asymptotic  expansion  is 
the  likely  thing  to  expect. 

The  equations  for  the  zero   order  approximation  are 

'^or  ~     oo-^o^r   o^^oa- 

o  =  /f^.+p^,^        y  o<>i<i 

(1+.8) 

IT  =   p       p    ,  _p       p        "^ 

0  =  /P    ,.+P   ^  >  1  <  ^l  <  1+r 

Vo^=l 

and  since  the  vorticity  is  zero 

If  V7e  assume  u  (-oo,'i)  =  U  (-co  j-^i,)  =  1,  the  solution  of  this  set  is 


;■  >  C! 
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%  =  1 

'^0  =  1 

'■o^'^ 

F„  ='l 

p„  =  /(1+pr-*) 

?,  =  lf 

('-L.9) 


The  zero   approxlv.iat Ion  represents  a  parallel  flovj  vjith  velocity  c 
in  both  layers. 

The  equations  for  the  first  order  approxiraation  are 

0  =  l-/fl^-Pl/^  0  <  ^  <  1 

0  =  ^i.>i-^^i  y 

(i^.io) 

Uicr  =  -  '^'■^l^-Pl.r  ^) 

0  =  1-/Fi^-Pi>^  \  1  <  ^L  <   1+r 

subject  to  the  vortlcity  conditions 

^1^><.=  °  '       ^1>5.1=  °  • 

The  solution  of  this  set  of  linear  equations,  subject  to  the 
appropriate  boundarjr  conditions  from  (i)..&)  is 

f]_  =  &^{^)^.  F^  =   A^(o-)(^-l)+a^(a-) 

(ll.ll)        u^  =  -a^(<r)  U-^  =  -A^icr) 

p^  =  -/a^(o-)  +a^(.T-)  P^^  =  -jrA-L(<r)(/?.-l)-/a^(o-) 

+(a-l-pr)+k^  +A^(cr)  +  (\-l-r)+K^ 

where  k,  and  K^  are  constants;  and,  due  to  the  boundary  conditions 
referring  to  pressure  at  'c  =  1  and  ^i=   1+r,  a,  (r)  and  A,  (o^)  must 
satisfy 


-  .1 
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(l-/+,o/)a^(<r)-pA^(cr)  =  P%-k^ 

/a3_(oi+(/r-l)A3_(cr)  =  K3_ 
or 

(l-/+p/)ai(.-)-pAj^(cr)  =  0 

(1+.13)  ,     ,  , 

/a^(:T')  +  (/r-l)A^(o-)  =  0  . 

It  is  here  that  the  bifurcation  phenomenon  appears.   If  each  of 
a,  (cr)  and  kAo')    is  to  be  non-zero  then  /  must  satisfy 

(ll.ll;)  (l-/+p/)(/r-l)+p/  =  0 

or 

(Lt.l5)  r(l-p)/2-(l+r)/+l  =  0 

which  gives 


/  = 


"SFTI^pT 


These  values  of  /  are  just  equal  to  the  values  X,  =  gh/c,  ; 
Xp  =  gh/cp  which  determine  the  critical  speeds  of  the  linear  theory, 
as  found  in  the  preceding  section.   In  order  to  obtain  motions  other 
than  parallel  flow  vre  assume,  of  course,  that  /  satisfies  (Li.,l!|). 
The  functions  a,  (cr )  and  A,  (or)  cannot  be  determined  at  this 
stage;  to  find  them  vie  must  pass  to  the  set  of  equations  obtained  by 
carrying  the  development  to  next  order.   The  equations  for  the 
second  order  approximation  are 

^20-  =  a-j_(.x)'i-a3_(o-)a^(o-)-/f2^-P2,j-~^ 
([|..16)    a^(cr)4  =  a^(cr)-/f2/,-P2^;.  ('       0  <  ^<  1 

0  =  f2,,+  u^-a^Cr) 


,  "^  ""^  !t?  '*'■  ■'*'■  "*"       'V  •  •  r>  fV  "     I 
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(>)..16)        U2^=  A[{(^){^-l)+a[{cT)''AJ_i(r)A[{'J)-jtF^^'?^^\ 
A![((r)U-l)+a^(/r)  =  A^(<r)-/P2,j-P2^_  \      l<^^<l+r 

subject  to  the  vorticity  conditions 

and  the  boundary  conditions  from  (k.6).   The  solution  of  this  linear 

system  is 

3 
^2  =  a2(<r)/c-aT^(c-)  -^ 

2  "     ^ 

Up  =  a^(cr)-a2(vr)+a-|^(<r)  -^ 


P2cr  ~  ^2^^^"^'^l^^^^l^'^~^'^^^l^'^^"'^^2^'^^^'' 


(i^..l7) 

P2  =  -A^(a)il^  -a^((r)i^-^  +A2(o-)  (/?-l)+a2(o-)-  -\- 

U2  =  Aj(a-)-A2(cr)+a^(cr)('J-l)+A^(j)  ^'-p- 

P2^=  A2(a-)-/a2((r)-3A-L(cr)A^(a-)+a^(a-)+  ^  a'^'(o-) 
+  [A^(o-)-/A2((r)-a^'(--)](A(.-l) 

.  [/a>)-A';(^-)]  ^^}l  ^/A';(a)  i:^ 

provided  that  the  functions  ap(<r),  Ap(c-)  satisfy  the  ordinary  dif- 
ferential equations 

(l-/+p/)a'(<r)-pA'(<r)  =  I 
(i|.l8)  ,     2     ,2 

Xa2(.T)+(/r-l)A2(a-)  =  J  . 

The  equations  ([j.,18)  coiae,  of  course,  from  the  boundary  conditions  on 
the  pressure  at  the  free  surface  and  at  the  interface.   In  these 
equations  the  I  and  J  are  given  by 


(  ^s'^' 


1-M>, 
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I  =  ^(3-/+p/)a"0r)  +  (p-i)a|(cr) 
+  3a^(cr1a^(a-)-3pA^(o-)A^((r) 

and 

J  =  |(/-6r+3/r2)a'^(cr)+^(/r^-3r2)A"^(cr-) 
+  a^(:r)+rA^(o-)-3A^(T)A^(cr)   . 

Thus  I  and  J  are  defined  b'/  the  first  order  terms.  \Ie   observe  also 
that  A'((r),  for  example,  is  determined  from  (Lj..l3)  as  soon  as  a^(o-) 
is  known.   Since  the  deterr-^inant  of  (i|..l8)  is  zero,  as  we  know, 
we  have 

(!;.19)  (/r-l)I+pJ  =  0 

as  an  equation  for  the  determination  of  a-,(c5")  and  A^((J').   If  we  use 
the  relations  (!.[.. 12)  and  ([f-.lli.)  the  equation  (i|..ll4.)  reduces  to 

(I(.,20)  mQa^(o-)  =  m-|_a^(5-)a^(a- )+m2a^  (cr) 

where 

m^  =  ^  [r(r+l)/-(r2+l+3pr)] 

rn^  =  Ya^   [/2(r2.1)-h/(l-2r)+l] 

m^  =  -[2/r(p-l)+r+l]+  ^^  \    . 

The  general  solution  of  ('i,20)  is  expressible  in  terms  of  the 
Jacobian  elliptic  functions  usually  denoted  by  en  and  sn  ;  and 
these  lead  to  cnoidal  v;av6S.  For  these  waves  the  quantities  h  and  H 
must  be  Interpreted  as  mean  depths  of  the  layers  and  c  must  be 
interpreted  as  a  mean  velocity.   However,  in  what  follows  we  vjill 


o;\!(i 


o:. 


O'l 
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confine  our  attention  to  the  solitar^.^  x-jave.   For  the  solitary  wave 
special  boundarj'-  conditions  raust  be  litiiDOsed.   Since  the  raediura  is 
supposed  to  be  at  rest  at  infinity  the  horizontal  velocity  there  vrith 
respect  to  the  movinp  coordinate  system  is  c  which  means  that  for  the 
dimensionless  velocities,  u(-oo,'t)  =  1  and  U(-co,'^)  =  1.   The 
velocities  u  and  U  satisfy  these  conditions  and  hence  i-ie   can  take 

0        0 

u,  (-oo,'2)  =  U,  (-co  ,/?)  --  0  which  raeans  a,  (-co)  =  A^(-oo)  =  0,  From 
{l\.,12)   we  see  that  the  horizontal  velocity  conditions  imply 

k^  =  0  ;  K^  =  0  . 

The  vertical  velocity  in  the  lower  layer  is  given  by  (ij..7)  and  using 
what  we  have  found  for  f  and  f, 

O        X 

V  =  e  y^  a^(Q")'l+( 


•  •  •     • 


The  condition  that  the  vertical  velocity  is  zero  at  -co  is  satisfied 

in  the  lowest  order  approxlnation  to  v  by  taking  a-j^(-co  )  =  0.   Since 

the  vertical  velocity  at  the  crest  of  the  wave  is  zero  we  take 

a  (0)  =  0,   In  addition  to  the  above  conditions  on  a,(a")  we  satisfy 

ti 
the  condition  u^(-co,^)  =  0  by  taking  a2(-oo  )  =  a^(-co  )  =  0  as  we 

can  see  from  ([{..17).   Subject  to  these  conditions,  the  integration 

of  ('4.20)  gives 


n      rn-,  a^  {-r) 


(1J..21)  m  a^(cj)  =  ~^^ +rapa^  (cr)  , 


If  we  multiply  (!^.21)  by  a^(/r)  and  integrate  we  find 
and  from  this 


C  '■ 


nriJ    -io 


r.-,l 


f&noo    "■ 


•eo.M 


^'. 


0   = 


ai    Of' 
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1  O 

The  function  AAc")    is  then  given  by 

If  we  now  suppose  that  what  we  have  found  gives  a  sufficiently 
accurate  representation  of  the  solitary  wave;  if  we  use  the  symbol 
='^  to  denote  approxiraation;  and  if  we  return  to  the  x,y  variables, 
a  collection  of  the  results  ,!?ives  the  following: 


r(x,Y)  =  h^t-  — -  (/-X)h>L  sech^  ^J  ^   {J(-\) 

u(x,y)  =  0+  --2  (/-A)c  sech^  "zE^rr  ^'^'"^^ 
1  ^    '  o 

o  3ni^  o  ^   /M^ 

p{x,Y)  =  5c2[x(l+pr-/^)-  ^  (/-X)(l-/'^)  sech'^  |e  ^  E^  (/"^^l 

v(x,Y)  =^c(/..A)2  ~2  /.^^  sech^  JL  y  _2(/.x)  •  tanh  ^7  ^(/-X) 

1    o  o  o 

3M,(/-X)h  ,      p  ,   /ra   , 

!J(x,Y)  =■  =-  ^  11^  /c  =ech2  |jy1i(/-U 


In  these  results 

c 
where  c  is  the  velocity  of  the  solitary  wave;  6  and  /\   are  the 
respective  densities  of  the  lower  and  upper  layers;  h  and  H  are  the 


■ion  0 


:\-: 


.Yi>-)« 


31j. 


resoective  depths  of  the  lovrer  end  uooer  layers  at  infinity.   In 

terms  of  the  vertical  distance,  y  ,  frova  the  bottom  at  x  =  -oo  , 

'  ■-  o 

Zo  .i0 

Y  =  j   t  (-oo,y)dy  =  J  cdy  =  cy^  . 
b  0 


The  constants  m.  are  {^iven  by 


=  1  r-nr-n-n  )   /^-C-n^. 

o 

1  "  TT?^  L/2(r2.1)+/(1.2r)+l] 


m_  =  -J  [r(r+l)/-(r'^+l+3pr)] 


m 

ra^  =  -  [2Xr(p-l)+l+r] 


and  /  is  a  zero  of 


r(l-p)/2-(l+r)/+i  =  0 


that  is 


(/-,/-   (l+r)+y(l-r)^+l4.rp 


or 


l  -    It     _  1+r-  y(l-r  )  '+l!.rp 
^        ^2  ^rd-pl      • 

It  will  be  noticed  that  there  must  be  a  resbriction  on  /-X  since  its 

m  (/-X) 

sign  must  be  such  that  -~- >  0.   This  quantity  appears  under  a 

o 
square  root  sign  in  various  places  above.   It  can  be  shown,  however, 

that  i''ip/ra  is  positive  so  th£.t  the  condition  becomes  /-X  >  0. 


5.   D is cuss  ion  of  the  Results 

According  to  the  approxiirjate  expression  for  f(x,Y)j  the  maximum 
elevation  or  depression  of  the  solitary  wave  at  the  interface  of  the 
two  layers  is 
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q   =    !r(0,ch)-h{ 
and   frci  this 
(5.1)  +q   =    -  -~^    (/-X)h    . 

In   terr.s   of   this   a-'aplitude   the   velocity  of  the   wsve    Is    given  by 

.qh  _  ,   0^   ,/       q    ^'-1 
2  -   ^  -  ^   ±hW; 


c 
or 


^^'^'  ^  -  /  1  q  hUr-l)T2/r(p-l)+l+rT  * 

The  strean  line  functions  can  be  ex-ore ssed  in  the  form 


f  (x,y)  =  hi  +  q^   sech^  ^^y+q(_l.) 

o 

(5.3) 


where 


p<.,„  =  M  i ,  liii^i^^  .och^  iyi,(^^) 


"^L  =  3[r(i-^-l)+/(l-2r)+l] 


'-'%        (/r-1)  [r(r+l)/-(r^+l+3pr)] 
and  the  sign  associated  vjith  q  is  to  be  chosen  so  that 

-^  o 
It  can  be  shown  that 

(/r-l)[r(r+l)/-(r^+l+3pr)]  >  0  . 

Hence  the  si.n-n  to  be  taken  is 


^Krt     -.r 
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and 


+  if  /^(r^-l)+/(l-2r)+l  >  0 


-  if  /^(r^-l)+X'(l-2r)+l  <  0  . 


The  possibility  of  having  /  (r  -l)+/(l-2r)H-l  =  0  will  be  considered 
later. 

The  solitary  wave  for  a  single  layer  of  depth  h  and  density 
equal  to  the  constant  5  can  be  obtained  from  the  above  by  setting 
r  =  0.   This  gives 


/  =  1  , 


-3ra, 


3  >  0 


so  that 


r(x,Y)  =   h1+   q^  sech^  Ie^/Ii 

n(x,Y)  =    c    -    qgsech2|./5 

p(x,y)  ^5Eh(l-/0  +  5c2  g   (1-/-/)    sech^  J^/S 

vIX,y;  Cs/^     ,^    oech     211'^^     ^anh  -^  y  ^ 


and  the  velocity  of  the  wave  from  (5.1)  is 


"h  /^  ,    q 


or 


gn       H 


which  is  the  Scott  iiussell  formula.   The  stream  lines  are  given  by 


y  =  f(x,Y)  =  ^3^  =  ,73  cfe  =  TT  ^ 


where 


37 
iS.k)  ^3  =  h  +  q  sech^  1-,/^ 

is  the  equation  of  the  profile  of  the  solitary  x-rave  at  the  free  sur- 
face.  In  terms  of  y  the  horizontal  velocity  corarjonent  is 

(5.5)  S(x,y)  =  u(x,y)  =1  (2h-yg) 
and  the  pressure  is 

p(x,y)  =  p(x,y) 

=  5g[h(l-1)  +  (l+g)q(l-'?)  sech^-l^-V^  ] 

(5.6)  p'(x,y)  =  5g(yg-y) 

which  shovjs  that  the  hydrostatic  lax^j  holds  up  to  second  order  terms 
in  q.   The  approximation  to  the  vertical  velocity  component  is  of 
higher  order  than  the  approximations  to  the  other  quantities.   It  is 

v(x,y)  =  v(x,y) 

-  -cy-j^  —-2  (yg-hJ  tanh  ^y-|^   . 

The  results  (3,lj.),  [y-^),    (5.6)  and  (5.7)  coincide  with  those  given 
by  Keller  [151. 

For  the  case  of  the  two-layer  medium,  /  may  take  either  the 
value 


"  ^2  ~  2rTl^^pl 


l+r+y(l-r)'^+lfpr 
which  is  positive,    or   the   value 


V  -   i/     =  l+r+../(l-r)^+l|.pr  ^ 2 

^        '^l  2rTr^^     ■  ; 

l+r-y(l-r)'"+4pr 


."■:)c  - 


(  . 


"^?5T0    h'-'no" 


rfi?r:i    z\ 


il 


lo'iqqa   sa'i.'      .p  ni 


:t.c; 


2::? 


;f<''T' 


•  =      ..X 
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which  is  also  positive  because  p  <  1,   These  values  lead  to  two  very 
different  kinds  of  solitary  wave  behavior. 
Type  1.   If  /  =  /2>  ^^en 


^1     _  3[(/2r-l)V2(l-/2^^ 


and  since 


"^"^o        a-j(^r)[3pT-^il'rJ(^)+r^{l"J(^)] 


1_/  r   =  l-r+v/(l-r)'^+[{.pr     ^      q 
1  +r  +  y  (l-rT^+lf.  pr 

1-/       =  -(l-r)+v/a-r)^+Ij-pr     ^   ^ 
l+r+  y(  1 -r  )  ^+[|.pr 


•^ —  is  positive  and  the  sign  to  be  associated  with  q  in  (5.1),  (5.2) 

~  '  o 

and  (5.3)  is  the  positive  sign.   If  xje  introduce 

(/pr-l)2  +  /.(l-/p) 
(J   =  £ ^ S >  0 

^   3pr+(l-r/2)+r2(l-/2) 
we  see  frora  (5.3)  that  the  strean  lines  are  given  by 


.2  X     /    ^^^1 


y  =  hn+q^c  sech  ^7^(1./  y  0  <  ^t  <  1 


and 


A,' 


[/^(^.-l-r)+l]     2  X  /  ^"^^-1 
J   =  h'i+q  — -i^TT^p-^ sech  ^^y-.^-^-^     1  <  /^  <  1+r 

and  we  see  that  these  are  lines  of  elevation  above  the  stream  lines 
of  the  uniform  flow.   The  wave  profile  at  the  free  surface  is  given  by 


p  y   /  3qd-, 
y^  =  h(l+r)+  -^-|-j,  sech  ^J^^^ 


2 
or 


y^  =  h(l+r)+Q  sech"^  "SI  ^^ 
where 
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Q  =  rr/;F 


2 

r 


is  the  araplitude  of  the  wave  at  the  free  surface.   The  wave  profile 
at  the  interface  is  given  by 


y^  =  h+d-Z^rjQ  sech^  -g^  7  -^i  . 

The  amplitude  of  the  wave  at  the  free  surface  Is  l/(l-/pr)  tii-nes  the 
amplitude  at  the  interface.   If  the  densities  of  the  layers  are  nearly 
equal  this  factor 


1   _  l+r+  y(l-r  )'^+ij.pr 
''    l-r+ /(l-r)'^+ll.pr 

is  approximately  equal  to  1+r.   If  the  depths  of  the  layers  are  the 
same,  the  amplitude  at  the  free  surface  is  about  tv.'ice  the  amplitude 
of  the  wave  at  the  interface. 

The  wave  velocity  when  /  =  /p  is  given  by 

gh  ^y  ,  ,  Q   (/gr-l)^fp(l-/,) 
-^  -  h"-  "E   UXprd-pj-d+rjJ    • 

For  the  case  of  layers  of  equal  depth 

r  =  1         /p  =  — 1- 

1+yp 

and  the   ijave   velocity  is   given  by 


(  r.-^-  r' 


Vi  i."ja   Q'i.''   C:i;:73i 


'xo^too'l   c.cfiJ   I 


•^    -^'i  r'  <■    ■" 


ko 


£h  *^  __1 Q_   .   1__ 

This   coincides   with  the   Scott  Russell   fonnula   for   p  =  1. 
Type   2.      If  /'  =  /^  the   sign  of 

/^(r2-l)+/^(l-2r)+l   =   (/^r.l)2./^(/^-l) 


=  ■4   [r^-r+2pr+(r-2)  y(l-r)^+lj.pr] 


may  be  positive  or  negative.   It  is  negative  for  example  if 
0  <  r  <  1;  anci  it  is  positive  if  2  <  r.   In  either  case  the  amplitude 
of  the  wave  at  the  free  surface  is  l/|l-/,r|  times  the  amplitude  of 
the  wave  at  the  interface.   This  factor  has  the  value 


0  <   1   =  r+l-y(l-r)^+ltpr  ^ l,ij(l-n) 

■^  r-l+y(l-r)^+!ipr   (1+r)  [r-l+ y'Cl-r  )2+!;pr] 

and  we  see  that  if  the  densities  of  the  layers  are  nearly  equal,  the 
factor  is  small.   If  r  =  1  the  factor  has  the  value 

1     _  l-s/p 

The  waves  of  the  present  type  tend  to  have  quite  small  amplitudes 
at  the  free  surface.   They  are  internal  solitary  x^faves,  and,  as  we 
will  seem  a  moment,  the  vave  at  the  interface  may  be  a  wave  either 
of  elevation  above  the  undisturbed  line  or  of  depression  below  it. 
In  connection  with  the  above  remark  about  the  amolitude  at  the  free 
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surface  it  should  be  noted  that  there  is  a  strea^n  line  in  the  upper 
layer  which  is  -parallel  to  the  bottom.   As  vie  can  see  frora  (3.3) 
the  horizontal  streara  line  occurs  when 

1  =  1+r  -  -^ 

=  l+r+  yjl+r  )^-J-|-(l-p)r 
2 

When  the  density  difference  is  small  this  line  is  near  the  free  sur- 
face . 

If  /^(r  -!)+/,( l-2r)+l  is  positive  the  positive  sign  must  be 
taken  in  (5.3)  and  the  strea^n  lines  are 


y  =  h',+q>i.  sech^  ^Jc^{-^^) 


^-o' 


0  <  'L  <  1 


and 


y  =  h>L-q     X'^^.i sech'^  ^yq(-^)    i  <  /^  <  1+r, 

These  are  lines  of  elevation  for  0  <  A  <  1+r  -  -r-  .   For 

1+r  -  J-  <  ^l  <   1+r  the  stream  lines  in  the  upper  layer  are  lines  of 

denression.   The  internal  wave  of  elevation  at  the  interface  is 


^         2  X  /    "^1 

yi  =  h+q  ^ech  ^^yq(-^_^) 


and  the  small  amolltude  wave  of  depression  at  the  free  surface  is 


J,   =   h(l+r)  -  j^~^   sech^  |^yq(33^-^ 


)  . 


k2 


If  /^(r   -1)+/,  (l-2r)+l   is   negative   the  negative    sign  raust  be 


taken   in   (5.3)    and   the    streava  lines 


y  ^  h'l-q^   sech^  _^^yq(-j^.l^)  0  <  ^^<  1 


o 
and 


y  =  h»t+q  — L-p— ^^ sech''  ^/q(:p^)  1   <^<  1+r. 

These  are  lines  of  depression  for  0  <^  <  1+r  -  j—  .  In  the  remaining 
part  of  the  upper  layer  the  stream  lines  are  lines  of  elevation.  The 
internal  solitary  ucve  of  depression  at  the  interface  is 


2  X  /   ™1 
y.  =  h-q  sech  ^yq(^^) 

and  the  small  amplitude  wave  of  elevation  at  the  free  surface  is 


^s   =   h(l^r)+q(-^^)  sech2  ^V^i^^)    . 
For  a  medium  in  which  the  depths  of  the  layers  are  the  same 


r  =  1  ;         /-,  = 


'^    1-yp 

and 

/^(r2-l)+/3_(l-2r)+l  =  1-J(\   =  -=-^^  <   0  . 

1-yp 

Therefore  in  this  case  tho  internal  solitary  i\Tave  at  the  interface 
is  a  depression  wave. 

The  wave  speed  for  x-javes  of  Type  2  is 
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air^BP. 
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72  -  '^li^  h(T{F-TTT2:f]T(  p-i )  +1TFT  * 

Since  /]_i"-l  >  0  and  P/,r ( p-1  )+l+r  <  0,  the  formula  for  the  speed  can 
be  written 

.h        ,  !/^(r2-l)+/^(l-2r)+l| 


■2  -  /i-q  Errji^Tr^Tj^TT^pT-Tr+FTT     • 


For  r   =   1   this    is 


gh 


2-:^ 


c   "  1-  yp     2h  yp 

1-  /p  2h  /p 

or  to  the  same  order  of  approximation 

Let  us  suppose  now  that 


//  _  /•  -  l+r+y(l-r)'-+'4pr  _ 

A  -  Ax  rFir^pT 


l+r-.y(l-r  )^+lj.pr 
and  that  the  depth  ratio  r  and  the  density  ratio  p  are  such  that 

f  , 

/^(r2.1)+/^(l-2r)+l  =  -^  [r2-r+2pr+(r-2)7(l-r)^+l|pr  ] 

vanishes.   In  order  to  discuss  this  situation,  for  which  m,  is  zero, 
V7e  need  to  return  to  the  different ic.l  equation  which  determines 
a-(o").   If  m,  =  0,  the  equation  is 


kh 


'     2      2 


The  general  solution  of  this  equation,  subject  to  our  conditions  at 
infinity.  Is 


I     ^  ra 

lae  cr  >   0 


^i'^>  =  \  .3 


Pe   °  (T  <   0 


\ 


This  solution  is  continuous  with  a  continuous  first  derivative  if 
and  only  if  each  of  the  constants  a  and  p  is  zero.   The  choice  of  the 
latter  value  for  a  and  p  leads  to  f^  =  0  and  we  are  left  with  a 
parallel  flow. 

It  mi£-hG  be  supposed  that  the  reraaining  solutions  may  lead  to 
something  more  interesting.  For  example,  the  discontinuous  solution 

a,  C;r)  =  aixicr)e  +   Pu.(->J-)e 


(where  a  f^  p  and  [i(cr)  is  the  unit  function)  suggests  that  internal 
bores  may  exist  at  the  interface.   If  we  apply  the  first  of  the 
shock  conditions,  which  requires  that  mass  be  conserved  across  the 
discontinuity,  v/e  find  (using  the  approximations  under  consideration) 
that  p  must  be  equal  to  the  negative  of  a.   Thus  the  first  shock 
condition  requires 


a-,  (o")   =  a(si2n'xii  a  ■)    e 


ill 


julor 


1+5 

Now  if  ve  apply  the  second  shock  condition,  vjhich  requires  conserva- 
tion of  momentun  across  the  discontinuity,  it  turns  out  that  an 
additional  relation  between  r  and  p  riust  be  satisfied.   This  relation, 
however,  is  incompatible  with 


'^-r+2pr+(r-2)  Jil-vW+hrp   =  0  , 


and  hence  a  bore  with  the  ^^roperties  our  theory  requires  cannot  be 
admitted. 


The  only  remaining  possibility  is  to  take 

-kiy 


a^C^r)  =  ae 


J2 

ra 


but   this  gives   a   solution  for  which  the  vorticity   is   not   zero. 
It    seeras   then  that   for   the    special   circurastance 


'^-r+2pr  +  (r-2)  y(l-r)'^+i|rp  =  0 


the  only  flow  of  the  second  type  vjhich  our  theory  allov:3  is  one  of 
parallel  flow. 
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Part  II 
Liquids  ^^7ith  a  Continuous  Vai^iation  in  Density 

1 .   Introduction 

This  part  presents  an  analysis  of  the  problera  of  solitary  waves 
in  an  incompressible  mediura  ^^Jith  a  continuous  variation  of  density. 
The  general  theory  is  developed  under  the  assuroption  that  the  density 
variation  is  arbitrary,  but  in  solvinr'  concrete  probler.is  we  assume 
that  the  density  of  the  nedium  in  the  equilibrium  state  increases 
exponentially  vjith  the  depth  of  the  mediura  in  the  channel.   The 
method  we  use  is  applicable  to  other  kinds  of  variation  of  density, 
but  these  in  general  would  lead  to  the  necessity  of  solving  certain 
ordinary  second  order  linear  differential  equations  with  variable 
coefficients  instead  of  those  with  constant  coefficients  which  arise 
for  the  case  of  exponential  variation  of  density. 

Although  the  analysis  in  this  part  is  similar  to  the  analysis  of 
Part  I  for  solitary  waves  in  a  medium  havinri  tvro  layers,  there  is  a 
mini'^-um  of  reference  to  Part  I  and  hence  the  development  of  this  part 
is  virtually  self  contained.  The  general  nonlinear  problem  is  formu- 
lated in  Section  2.  In  Section  3  the  nonlinear  problem  is  linearized 
and  it  is  shown  that  the  linear  theory  furnishes  an  infinite  number 
of  critical  speeds. 

In  Section  l\.   there  is  a  return  to  the  nonlinear  problem  and  it: 
is  analyzed  by  using  a  variant  of  the  shallow  v;ater  technique  invented 
hj   Priedrichs.   The  analysis  shows  that  solitary  wave  solutions  (and 
cnoidal  wave  solutions)  appear  as  bifurcations  from  states  of  uniform 
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parallel  flow  with  speeds  coinciding  with  the  critical  speeds  found 

in  Section  3.   The  solitary  wave  solutions  can  be  denoted  by  S 

where  the  subscript  n  indicates  that  the  solitary  v;ave  defined  by  S 

travels  with  a  speed  in  the  neighborhood  of  c  ,  the  n   critical 

speed  of  the  linear  theory.   The  critical  speed  c  approaches  zero 

as  n  — >  00  ,  and  in  fact  the  set  J  c  v  is  a  discrete  set  of  mono- 

(  nj 

tonically  decreasing  values  which  has  zero  as  its  only  linit  point. 
The  solitary  vjave  derived  frora  S,  travels  with  a  speed  in  the  neigh- 
borhood of  the  highest  critical  speed,  and  it  hes  'oroperties  similar 
to  those  of  the  ordinary  solitary  wave  in  a  medium  of  constant 
density.   The  solutions  S   (n  7=  1)  are  found  to  yield  internal 
solitary  waves.   The  amplitude  of  each  of  these  waves  at  the  free 
surface  is  small  compared  with  the  maximum  departure  of  certain 
internal  stream  lines  from  the  parallel  lines  ii/hich  they  approach 
at  infinity.   The  state  of  the  medium  for  these  waves  can  be  described 
roughly  by  saj^ing  that,  starting  at  the  bottom,  the  channel  is  divided 
into  adjacent  horizontal  strips  of  nearly  equal  widths  (except  the 
uppermost  one)  vjhich  alternately  contain  only  stream  lines  of  eleva- 
tion and  only  stream  lines  of  depression  with  respect  to  the  undis- 
turbed parallel  flow,  of  course.   The  horizontal  division  lines  for 
the  strips  are  themselves  stream  lines;  the  highest  such  division 
line  is  found  to  lie  a  small  distance  below  the  free  surface;  and 
there  are  n  strips.   The  maximum  deoarture  of  the  stream  lines  from 
the  equilibrium  lines  of  parallel  flow  in  any  one  of  these  strips 
(except  the  one  which  contains  the  free  surface)  is  slightly  greater 
than  the  raaximum  departure  of  the  stream  lines  in  the  strip  immedi- 
ately below.   These  proioerties  and  others  are  discussed  in  some  detail 
in  Section  5,  where  particular  attention  is  paid  to  the  waves  derived 
from  S,  and  Sp, 
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2.   Formulation  of  the  Problerg  for  the  Case  of  Variable  Density 

For  the  purpose  of  studying  tvro-diraensional  waves  of  permanent 
type  in  an  inviscid,  incorapressible  rr^edium  of  variable  density  vje 
suppose  that  a  cross  section  of  the  medium  in  the  equilibrium  state 
fills  an  infinite  horizontal  strip.   The  medium  is  assumed  to  be  sup- 
ported by  a  rin'id  bottom,  and  the  uoner  surface  is  assumed  to  be  a 
free  surface  at  which  the  pressure  is  zero  end  there  are  no  geometric 
constraints.   We  also  suppose  that  a  disturbance  of  some  kind  in  the 
medium  initially  at  rest  has  created  a  wave  of  permanent  type  which 
moves  to  the  left  with  constant  velocity  c;  and  that  the  character 
of  this  wave  is  observed  from  a  coordinate  system  which  moves  with 
the  wave.   VJe  take  the  x-axis  of  this  system  so  that  it  coincides 
vjith  the  bottom  of  the  channel  and  the  y-axis  so  that  it  is  positive 
upward  and  passes  through  a  crest  of  the  wave.   V'ith  respect  to  these 
axes  the  wave  is  stationary;  and  the  velocity  of  the  medium  at 
infinity  is  horizontal  and  equal  to  c.   V'e  take  the  depth  of  the 
medium  at  infinity  to  be  h. 

The  differential  equations  for  the  flow  under  these  conditions 
are  well  knoiNm  to  be 

5(uux+vuy)  =  -p^  , 


(2.1) 


6(uv^+vVy)  =  -6s-p^^  , 
?  -  Vx"'^v  ' 


u  +v   =  0  , 

X  y  ' 

us^^+vsy  =  0  , 

u5^^-Hv6y  =  0  , 


liJ. 


where  5(x,y)  is  the  density;  u(x,y),  'v(::,y)  are  the  horizontal  and 
vertical  velocity  components;  p(r.,y)  Is  the  pressure;  g  is  the 
acceleration  due  to  gravity;  ^(x,y)  is  the  vorticity;  and  s(x,y)  =  0 
is  the  equation  of  a  boundary  curve. 

If  the  stream  function  is  \l/(x,y)  so  that 

u^=  ^1/  (x,y)  ;  V  =  -il/^,(x,y)  , 

the  totality  of  stream  lines  is  given  by 

\!/(x,y)  =  Y 

where  y  ranges  over  a  finite  interval.   The  function  \!/(x,y)  can  be 
chosen  so  that  for  the  bottom 

Y  =  '!/(x,0)  =  0 

and  then  the  value  of  y  ^-'hxch  determines  the  free  surface  can  be 
taken  to  be 

Y  =  ch  . 

It  Is  assujaed  that  to  each  value  of  y  there  corresponds  one  and  only 
one  stream  line  which  is  given  by  solving  y  -   'l'(x»y)  for  y,  namely 

y  =  f (x,y)  . 

Let  us  use  x  and  y  as  independent  variables.   As  the  basic 
dependent  variables  we  choose  the  stream  line  functions  f(x,Y)»  the 
horizontal  velocity  component  u(x,y)  =  urx,f(x,Y)];  and  the  pressure 
P(x,y)  =  p[x,r(x,Y)].   ^^ith  these  variables 
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''    ^y   u  ^    ^y    u 

and  xije  find,  as  we  did  in  Section  2  of  Part  I,  that  the  first  three 
of  equations  (2.1)  reduce  to 

6u   =  -f  p  +f  p 

(2.3)  ^('^^xx-^Vx^  =  -^"sVPy 

!)   ^x     1   -^   ^■'"^x 
^       <?x   ^-  '        2  ^T   ^2  -"   • 
Y  Y 

The  fourth  and  fifth  of  equations  (2,1)  are  automatically  satisfied 
by  the  new  variables;  and  the  sixth  equation,  which  arises  from  the 
assumption  of  incompressibilitj'-,  becomes 


u5y+v5^^  =  u5^  =  0 


which  shows  that  0  =  5[x,r(x,Y)J  is  a  function  of  y  only.   VJith  the 
equations  (2,3)  we  see  from  (2,2)  that  we  also  have 

(2.i;)  af  Y  =  1  , 

In  order  to  make  the  equations  diriiensionless  i-ie    introduce  the 
variables 


5  =  ^ 

v=? 

?  =  |? 

'^=^ 

p  =  --a^ 

^        6 

S^c 

o 

0 

-  =  i 

^-'l 

c 
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where  5   is  the  density  of  the  raedium  at  the  bottom  of  the  channel, 
o 

The  equations  in  dimensionless  form  are 


6u^  =  -f/^P^+f^P,^ 
5(uf^^+u^f^)  =  -X5f^-p^ 
uf^  =  1 


(2.5) 


2 


where   5,    as  we  know,    is   a  function  of  /^   only.      These   equations  must 
hold  for   the    strip 

-co    <5<oo;  0</^<l» 

The  boundary  condition  at  the  bottom  is 

(2.6)  f(?,0)  =  0 

and  the  boundary  condition  at  the  free  surface  is 

(2.7)  P(C,1)  =  0  , 

The  equation  which  expresses  Bernoulli's  law  is  found  by  elimi- 
nating u  from  the  first  tliree  of  equations  (2,5)  and  then  integrating 
the  resulting  equation.   The  elimination  of  u  gives 

and  integration  along  ^l   =  cons' t.  gives 
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2 

1+f  r 

(2.8)  I  ( — ^)    +X5f+p  =  C(^)  . 

If  we  differentiate  the  Bernoulli  equation  (2.8)  with  respect  to 
/l  and  use  the  second  and  third  equation  of  (2.5)  we  find 

P  2 

1+f  f       1+f 

(2.9)  S  -I-  (—pi)  -  5  ^  (7i)+ot  [ 1  +  Af  ]  =  C  (>.)  . 

2  <?>i    f  2      ^^   1^       2f^ 

This  is  the  equation  vjhich  -must  be  satisfied  by  the  stream  line 
function  f(^,7.).   The  equation  can  also  be  interpreted  as  a  relation 
between  the  vorticity,  pressure  and  density: 


3«   The  Linearized  Equations.   Critical  Speeds 

The  equations  (2,5)  and  the  boundary  conditions  (2,6)  and  (2.7) 

are  satisfied  by 

u  =  1 

f  =^ 

1 

p  =  X  f  5(A)d/L 

which  expresses  a  state  of  parallel  flox-j  x/ith  constant  velocity.   Let 
us  consider  a  v;ave  motion  v/hich  is  a  sriall  disturbance  superposed  on 
this  state  of  constant  floxj.   Thus  vie   vjrite 

u  =  1  +  u"  , 

f  =  ^c  +   f'""  , 
1 

p  =  X  f  5(/^)d^  +  p'""  , 
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with  u' ,  f"'',  p"  regarded  as  small  quantities.   Then,  if  we  substitute 
these  quantities  in  (?.5),  ('?.6)  and  (2.7)  end  assume  that  the 
disturbance  is  so  small  that  second  order  products  involving  the 
starred  quantities  can  be  nefrlected,  we  find  the  linear  equations 


6Ug.  =  -Pr-X5fg. 

5f|^  =  -X5f^-p^ 
u'Vff  =  0 


subject  to  the  boundary  conditions 

f'""(^,0)  =  0 
p'"'(5,l)  =  0   . 

Prom  these  equations  we  find  that  f (^,^)  raust  satisfy 


(3.1)  ^^^'ie^-^^ln'i^-^^'  ^^l^-'^^l^  =  ° 


subject  to  the  boundary  conditions 

(3.2)  f|'(C,0)  =  0  , 

(3.3)  ^IV^'l)  =  Xf|'(?,l)  . 

Prom  here  on  we  assui'ie  that  the  density  decreases  exponentially 
as  we  go  upward  from  the  bottom  of  the  channel  to  the  top.  In  fact, 
we  assume 
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5(^)   =  e 
The   equation   (3.1)   then  becoraes 


(3.1|)  ^l^z''^l\n-^-^^^li-'^^'p  =  0   • 


This  equation  and  the  bou.ndary  condition  (3.2)  are  satisfied  by 


f ^  =  A  cos  v{E,+B)e^^   sin -^i  s/2kX-k^'V'^ 


where  A,  B  and  v  are  arbitrary  constants  except  that  A  must  be  small. 
In  order  to  satisfy  the  boundary  condition  (3.3)  we  raust  have 

(3.5)  (k-X)  tan  ySkX-k^-v^  +s/5a-lc-v^  =  0  . 
If  we  set 

(3.6)  a^  =  JzkX-k^-v"^ 

(3.5)  can  be  written 

2ka 

(3.7)  tan  a  =  -^ — ^ — ^     , 

aj+v--k2 

This  equation  relates  the  v:ave  nuiaber  v  to  the  speed  of  the  wave  and 
the  exponential  variation  of  density.   If  the  vjave  number  v  (or  the 
period  T  =  2'k/v)    is  given  fcr  a  fixed  depth  and  fixed  k,  the  equation 
(3.7)  has  an  infinite  nuraber  of  roots  a  ,  and  since 


-r'i  b; 
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there  are  an  infinite  number  of  speeds  corresponding  to  a  prescribed 
wave  length. 

In  Part  I  we  found  that  the  critical  speeds  could  be  defined  in 
various  ways.   Here  we  deterraine  them  as  the  speeds  corresponding  to 
waves  vjhose  xjave  length  tends  to  infinity.  For  the  case  of  exponential 
variation  of  density  vie   easily  see  from  (3.7)  and  (3.6)  that  these 
speeds  are  given  by  the  roots  of 

(3.8)  tan  a  =  — ^ — ^ 

where 


a  =  ySa^-k^ 


n 
and 

1  -  gh  _  a^+k'^ 

\-  y  — 2ir   • 

°n 

The  equation  (3.8)  has  an  infinite  nuraber  of  real  positive  roots,  and 
no  pure  iraaginary  root  of  magnitude  greater  than  zero.   It  furnishes, 
in  other  x^jords,  a  discrete  infinite  spectrura  of  critical  speeds.   The 
equation  (3.8)  will  be  discussed  in  greater  detail  in  the  next  section. 
For  the  time  being,  it  suffices  to  reriark  that  our  experience  with  a 
medium  of  tx-jo  layers  leads  us  to  suspect  that  there  exist  solitary 
and  cnoidal  waves  with  speeds  in  the  neighborhoods  of  the  critical 
speeds  given  by  (3.8). 
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k»      Solitary  Waves  Near  the  Critical  Speeds 

We  return  now  to  the  nonlinear  equation  (2,5).   For  a  variation 
of  density  in  the  forin  of  an  exponential  decrease  going  upward: 

0  =  8     ,  k>0. 


these  equations  are 


e  '^''"u^  =  -f^p^+fpp^ 


(l|.l)  e""^''''(uf^^+u^f^)  =  .Ae"2^\-p^ 

uf^  =  l 

subject  to  the  boundary  conditions 


f(?,0)  =  0 
(i^..2) 

P(^,l)  =  0   . 

We  proceed  to  investigate  a  solitary  vjave  which  moves  with  a 
speed  such  that  X  =  gh/c   is  near  some  value  /,  which  is  to  be 
determined,  but  which  will  be  found  to  have  one  of  the  values  fixed 
by  giving  c  any  of  the  critical  values  just  described  above.   Vie  use 
the  same  general  method  that  T:as  used  in  Part  I.   That  is,  we  first 
write  the  equations  (l|.l)  in  the  form 

(i|.3)      e-2k'L(^f^^+^^f^)  =  (/_^)  e-2i^\-Ce-2k\-p^^ 

uf,.  =  1  . 


57 


Then  vje  set 

/-X  =  e 

and  introduce  the  new  independent  variable 

The  equations  (i]..3)  thus  becorae 


e"^^^'u_=  -f.TD  +f  p^ 
ik.k)  ee"2^''^(uf^^+u^^f^)  =  ee'^^'^-Ze'^^^T^^-p^^ 

uf,,  =  1   . 


These   equations  must   hold   for 

-co    <c-<ooj  0<^2.<1 

and  the  boundary  conditions   are 

f(cr-,0)   =  0 
p(c^,l)   =   0      . 

We  assume  now  that  all  quantities  in  the  new  variables  can  be 
expanded  in  integral  powers  of  e,  thc.t  is,  we  assume 

00   , 
.CO..   V 

k^    ^ 

00.  y, 

P  (cr,  ^c,  e )  =  1_  e  P,,  (-^ ,  *i ) 

k=0    ^^ 


■X)    >  — .1   >     00- 
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th 
where  the  k   order  approxinations  are  subject  to  the  boundary 

conditions 

fj^(r,C)  =  0 


As  a  consequence  of  the  assuraed  expansions  ([j.,6),  the  vertical  velo- 
city must  have  the  expansion 

-   uf 

^  =  ^  =  -c-  =  ^f^  =77  uC 

=  -^   fW-^-^^Vo-'^Vlr^^---  ^  • 

¥e  intend  now  to  solve  the  equations  which  the  lower  order  approxima- 
tions must  satisfy,  and  which  come  from  substituting  {I^.6)  in  {k»h-) 
and  equating  coefficients  of  like  powers  of  e.   The  boundary 
conditions  (Lj..?)  are  set  in  the  ssxie  way  by  substitution  of  (!|.6)  in 
(I;.. 5).   The  expectation  is  that  the  analysis  will  show  automatically 
that  /  must  correspond  to  one  of  the  infinite  number  of  critical 

speeds  defined  'oj   the  linear  theory  of  the  previous  section. 

th 
The  equations  for  the  zero   order  approximation  are 

e  "^-^'u   =  -f  ^t)  +f  p  , 
('+.8)  0  =  /e-2=%,«„,^ 


subject  to 


f^(c',0)  =  0 


o 


p^(cr,i)  =  0  . 
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This  is  a  nonlinear  system  end  it  is  not  easy  to  find  its  general 

solution.   We  can  proceed,  however,  without  knowing  the  general 

solution.   Since  we  are  interested  in  solitary  xvaves  which  move  in 

a  medium  initially  at  rest,  we  can  assume  that  the  dirtiensionless 

zero   order  approximation  (  with  respect  to  the  moving  axes)  must 

exoress  a  state  of  oarallel  flow  with  constant  velocity  u^  =  1.   The 

,/       o 

solution  v;hich  expresses   this   state    is  readily  seen  to  be 


=      Ai 
o  <" 


u     =  1 
o 


_    /     ,    -2k'L     -2k  X 


VJith  this  zero   order  approximation  the  equations  for  the  first 


order  approximation  are 


e'^^a^  =  .^e'^^f^^ 


(4.9)  0  =  e-2I^^'^-/e-2k^f^^^.p^^ 

f^,^+u^  =  0 


subject  to 

f^(rr,0)  =  0 
P^(o-,l)  =  0  . 

This  is  a  linear  system  vjhich  can  be  solved  with  little  effort.   The 
elimination  of  u,  and  p,  fro:.i  (i;»9)  gives 

(i;.lo)  f3_-^_,,-.2kf^^,+2k/f3_^  =  0  . 


.:>.cfi   bail  c' 
£.5  ■:ixi:i   r 

r.l   5 von 


.ft      .J 


Tjt'i 


.  .  r  r  . .  5  J-  ^  . 


•<»ftxvori   snJ    c:    joe-':3- 


e^riT      .1  = 


n  To 


Vo     i:.-3^ 


■;[S^         '^71?  - 


•r/oe   srlT  noxo. 
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The  first  order  approxi-mation,  f-,  ,  to  the  stream  line  function  must 
satisfy  this  equation  in  the  strip  -co  <:j'<oo;  0<>t<l;  and  the 
boundary  conditions 

(^.12)  ^Ir^J'^'l^  =  /^l^C^,!)  . 

The  solution  of  (1|.10)  which  satisfies  (l|..ll)  is 


(i)..13)  f^,^=  a^('r)e^^'- sin  a^ 


where   a^  (cr)    is   arbitrary  and  a  j«^  0   is  defined  by 


(i^.ll^)  a  =  y2k/-k^      . 

(The  possibility  of  having  a  =  0  can  be  discarded  because  this  pos- 
sibility implies  a  state  of  constant  density,  the  analysis  of  which 
is  covered  in  Part  I,)  ?rom  (!i.l3)  integration  vjith  respect  to  '•^~ 
yields 

(I;.  15)  f^  =  a^(T)e^^  sin  a^L+h^ik) 

where  b^  Cv!)    is   an  arbitrary  function  of    ^L,      Since  vie   confine   our 
attention  to   solitary  waves,    x^e   impose   here   the   condition 

f^  (-co  /I)   =  0      . 


This  fixes  b-{^.).   In  fact,  it  permits  us  to  take  b,  (^)  =  0  since  v 


je 


choose   a,  (-co  )    =  0,      VJe  now  have 


61 

('4.16)  f^  =  a^('r)e^^sin  a>2.. 

[If  we  wanted  to  study  cnoidal  viaves  we  vjould  require  f,  to  be 

periodic  with  some  period  T-, ,  and  vre  would  also  require  that 
+T 
j   r^{r',^)dr    =  0.   This  xjould  fix  h^{i)    in  (!i.l5)  and  permt  us  to 


write 


^1  ~  3^1^^^)©^^^  sin  ah^ 


where  a,(cr+T- )  =  a,(a~)  and 


■-T+T, 
/■■  1 


a,  (cr')d'>  =0.]   Now,  in  order  to 


satisfy  the  remaining  boundary  condition  we  find  by  substituting 
(i|.l6)  in  ('-]-, 12)  that  we  must  have 

(i;.17)  a^(cr)[(k-/)  tan  a+a]   =  0  . 

It  is  at  this  point  that  the  bifurcation  phenomenon  appears  for  the 
case  of  variable  density.   We  find  a  solution  which  gives  more  than 
just  a  state  of  parallel  flow  if  we  take  a^{cr)   ^  0  and  choose  /  so 
that  it  satisfies 

(k-/)  tan  a+a  =  0 

or  what  is  the  same  thing 

(1;,18)  tan  a  =  —35 — ^     , 

a  -k^ 


■I    V' 


\      I    -^    »   T* 


n    ieu^ 


\  -a  »i- ; 
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The  bifurcation  values  of  '/ ,    in  other  x:ords,  the  values  of  /  which 
satisfy  (l|.l8),  are  just  the  values  \^   =  gh/c^  which  satisfy  (3.8) 
and  give  the  critical  speeds  obtained  from  the  linear  theory. 
The  equation  ('4..  18)  is  equivalent  to  the  two  equations 


tan 


a 


k 

a 


_  k 


-cot  -2  =   tan  (— ^)  =  - 


Hence  the  roots  of  (ij.,18)  are  given  by  the  Intersections  of  the 
hyperbola 


with  the  curves 


3  =  - 

p 

=  tan 

a 

p 

=  tan 

f^:^) 

These  curves  are  shoim  in  Pi!3.  'i.l.   Fe  see  from  this  figure  that 
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there  Is  a  denumerable  set  of  bifurcation  values  i^hich  increase  in 
magnitude.   If  k  is  snail  these  values  are  approximately  Integral 
multiples  of  -n:.   Stated  in  tei'ns  of  the  critical  speeds  of  the  linear 
theory  there  is  a  discrete  infinite  set  of  critical  speeds  which  has 
the  value  zero  for  its  only  limit  point.  We  suppose  hereafter,  of 
course,  that  /  is  a  bifurcation  value. 

The  first  order  approziination  to  the  horizontal  velocity  (from 
the  third  of  equations  (I}.. 9))  is  given  by 

^1  =  -  h'i 

('+.19)  ,. 

=  -a,('.-r)e    [k  sin  a-'/.+a  cos  a^_] 

and  the  first  order  approximation  to  the  pressure  (from  the  second 
of  equations  (m.,9))  is  obtained  from  an  integration  of 

from  which 

^  _  1  ,  -2k  ^-2k^v 
P3_  -  -^  (e   -e    ) 

(i^.20)  a,  (-T-)   ,        P 

~ —  e    [(a'^-k  )  sin  a^c   -21:a  cos  a  ^_]    . 

The  function  a,  (cr)  cannot  be  determined  at  this  stage;  to  find  it  we 
must  pass  to  the  set  of  equations  for  the  second  order  approximations, 
The  equations  for  the  second  order  approximation  are 

^"^^'  ^^27-  =  ^lrPir^2TPo/rPlT^l^-P2^r 
(4-. 21)  6    r^^-   e    fi>^-Xe    ^2^^2^i 

^2>f-'"lVL-'^2  =  0 


o:;3  3   erj.v  mo"'i) 


O-iJC 


subject  to  the  boundary  conditions 

f^ir^Q)    =   0 
P2(cr,l)  =  0  . 

By  elimination  of  u^  and  Pp  from  (!|«21)  and  use  of  the  first  order 
equations,  (!|,21)  leads  to 

(i;.22)  f2,^,,^,-2kf2^;^+2k/f2o-=  L(^'A) 

where 

+  3f;,;,(fl^^,.-2kf3_,^^)  +  2kf^^_ 

is   a  knovm  function  of  a,(a-);   and  fp   satisfies    (ij.,22)    in  the   strip 
-co    <o-<oo;    0<4<1  with  the  boundarjr  conditions 

(IJ..23)  f^{(r,0)    =  0  ; 

(,'4.21^)        f2T/?(^'^^   =  /f2,J'^,l)+3/fi^(r-,l)f-L^(cr,l)-f-Lj.(o-,l)    . 

The  equation  for  the  determination  of  a^  (o~)  cc.n  be  found  by 
multiplying  each  side  of  ([j..22)  by  e  "  sin  a'Z    and  integrating  the 
resulting  products  from  zero  to  one.   Th£.t  is,  if  v.e   integrate  the 
left  side  of 


1 

J 

0  0 


j  e'^''^  sin  a/i[f2,.^^^^-2kf2^2k/f2Jd^,  =  ]  e'^^  sin  a'l  L(^, /'Od>L 


00' 


C^.r) 


65 

by  parts;  use  the  boundary  conditions  (1^,23),    ('4.2!|)  and  assurae,  of 

course,  that  a  is  a  root  of  ('4.I8)  we  find  that  a,  (tr)  must  satisfy 

(1|.25)  -L 

e"^  sin  af3/f^/^(<r,l)f^^(o-,l)-f^,^(^,l)]    =  j   e"^'^'^  sin  a'lLi'T/Odt  . 

0 

Then  if  x^e  substitute  f,  =  a,  ((r')e   sin  a^    and  perform  the  integration 
on  the  right  hand  side  of  (!|.25)  it  turns  out  that  a^  (or)  raust  satisfy 

In  this  equation  the  quantities  m  ,  m. ,  and  nip  are  given  by: 

a   vk  +a  j 

(Z;.27)     m^  =   3/a(a  -t-k  )    [SCk^+a^je^^  cos  a-  (k^-a^)] 
■^   (k^-a2)(k"+9a2) 

m.  =  isflkfl  (a2+k2+2k)   >   0   . 
2    I^ka^ 

The  constant  m^   is  alvjays  greater  than  zero,  and  it  can  be  shown  that 


this  is  also  true  for  ra  .   Since 

o 


(k^-a^) 


cos  a  =  +  — p — p— 

k^+a'=^ 

the  constant  m^  may  be  greater  than  zero  or  less  than  zero  but  it 
cannot  be  zero  since  8e  +1  cannot  be  zero. 

The  equation  {L\.,26)   has  the  sai'ie  form  as  the  equation  (lj.,20) 
in  Part  I  which  determines  solitary  and  cnoldal  waves  for  a  medium 
of  two  layers,   Cnoldal  waves  for  our  medium  of  continuously  variable 


Y'i'^i.OisB     J.3>Jff:     ( 


:■•;.■;   n~.'V:  ^,fTi  bnB    ,  .n    , 
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density  can  be  derived  from  the  general  solution  of  (^.26).  Solitary 
waves  result  if  we  impose  special  boundary  conditions  on  ([i.26).  One 
such  condition  has  already  been  set,  namely  a.  (-co  )  =  0.  For  reasons 
explained  in  Part  I  we  also  impose  the  conditions  a-,  (-co)  =  0  and 

a. 
II 

a^(-co  )   =  0.      IJith  these    conditions,    an  integration  of    ('.i.26)    gives 

Tf  a-^  (cr ) 

(14-.  28)  m^a^(or-)   =  m^  — ^ —  +m^8i^{cr) 

and  another  integration  after  multiplication  of  (4.28)  by  a,  (?")  gives 

('4.29)  m^[a[{a-)]^  =  -^   aj(cr)  +m2a^(cr)  . 

Then,  the  solution  of  (Ij..29}  which  has  zero  slope  at  a- =  0  (conforming 
with  the  assumption  that  the  vertical  axis  passes  through  the  crest 
of  the  wave)  is 

3mp     p    /mT 
1  o 

¥e  suppose  now  that  what  vje  have  found  gives  a  sufficiently 
accurate  representation  of  solitary  waves  in  a  medium  with  the 
assumed  exponential  variation  in  density.   If  we  use  the  symbol  = 
to  denote  approximation;  and  if  we  return  to  the  x,y  variables,  a 
collection  of  the  results  shoiis  the  following: 


(/^.30)   f(x,Y)  =  h  -_£(/-x)he^^sin  a>L  sech^  ^J  ^ij(.x) 

1  o 


(1|..31)   u(x,y)  =  c-K^(/-\)ce^^(k  sin  a>i+a  cos  a1)sech^  "^^JT^'^'"^^ 

1  o 


I  Or     .: 


.:^) 
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(;4.32)   p(x,r)  '^  5^c2  f^(e-2^^-e-2k) 


v^ 


^  it^  e"^'''[(a2-k-)sin  a^-2ka  cos  a^]  • 


m^ 


(/+.33)  v(x,y)  =  c(/-X)2  --2  /-?  e^-^  sin  a't  • 


o  o 


In  these  results 


where  c  is  the  velocity  of  the  solitary  vjave;  and  k  >  0  is  the  index 
of  the  exponential  density  distribution  5{x,y)  =  5  e     .   In  terras 
of  the  vertical  distance,  y  ,  from  the  bottom  of  the  channel  at 
x  =  -00 , 

Y  Y 

lO  "^,0 


-  J      i|/y(-oo  ,y}dy  =   j      cdy  =   cy^ 


Th 


le    constants  m^  are   given  by   (1|..27)    and  /   is  to  be   determined  from 

2k/  =   a^+V? 

where  a  is  a  non  zero  positive  root  of 

.  2k  a 

tan  a  =      n   •■>:>      • 

a^-k"^ 

^2  m^ 

Since  — -(/-\)    aopears   under  the    square  root   siern  and   since  —  >   0   , 

o  o 

it  follows  that  /-X  must  be  positive,  which  means  that  the  speed  c 

of  the  wave  is  always  greater  than  the  critical  speed. 


s^ 


CO-    = 
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5»   Properties  cf  the  Solitary  V^aves. 
Solitary  'Javes  t-Jhose  Spe eds  are  n ear 
the  two  rlighest  Critical  Speeds . 

In  this  section  we  discuss  some  of  the  properties  of  the  solitary 
waves  which  are  obtained  from  (ii,30)  v;hen  k  is  small.   Particular 
attention  is  given  to  the  solitary  waves  vjhich  travel  with  speeds 
near  the  two  highest  critical  speeds. 

The  stream  lines  throughout  the  medium  are  given  by 


3m  •''m 

(5.1)  y  =  f(x,Y)  =  h>i--^(/-A)he^'^-  sin  a^i  sech^  ^/-^(/-X) 

1  o 

where  0  <  '^  <  1  ^nd  where  a  may  be  any   one  of  the  values  determined  by 

(5.2)  tan  1  =  1 
or 

(5.3)  tan(^)  =  -cot  I  =  I   . 


VJhen  k  is  small  the  values  of  a   which  satisfy  (5.2)  are  slightly 
larger  than  even  m.ultiples  cf  %,    and  the  values  of  a  which  satisfy 
(5*3)  are  slightly  larger  than  odd  multiples  of  it.   We  proceed  to 
discuss  these  tvjo  cases  separately. 

Case  I.   a  =  m2T[;  but  a  >  mSii  (m  an  integer). 

For  a  value  of  a  which  is  nearly  a  multiple  of  2%   we  have 

cos  a  =  ?5— -  =  — s — « 

l+tan"^  I   a'^+k'^ 

and  hence 


69 

^  =      (a^+-Ak)(a^-k^)(k^+9a^) 
^1    2a^(a^+k^)[8(k-+a^)e''^  cos  a-lk^-a^)] 

(a2+k^+k)(k^+9a^) 

= ., ^5 ^ =- 

2a^(a^+k'^)(ee^+l) 


which  is  positive.   If  we  vrite 


3rap 
<l=  -  ^   (/-X)  >  0  , 


^1 


the   streara  lines   are  given  by 

(5.i|.)      y  =  r(x,Y)   =  ht+qhe^^"^  sin  al  sech"  ^7  q(-J:) 

where 

,^  ^.  "^1    _         [Lka^(k^-i-a^)^(8e^+l)  

•^%        (k^+9a'^}[2k(k  -a^)+(k'^+a^)'^] 

In  terms  of  q,    a  and  k  the    speed  of   the  vjave    is   given  by 

m, 

or 


/c-  A1  gh  _  o^+k^n  iika^(8e'"^+l)q ^ 

c^  ^^^  (a^+k^+2k)(k'^+9a'^) 

Prom  (5.^)  we  see  that  the  stream  lines  are  parallel  to  the  horizontal 
bottom  of  the  channel  for  those  values  of  ^i   which  satisfy  sin  a^l    =  0, 
If  a  is  betv/een  zero  and  n  there  is  no  stream  line,  other  than  the 
bottom  one,  which  is  horizontal;  but  if  a  =  m2Ti:  where  m  is  a  positive 
integer,  the  horizontal  stream  linos  occur  xvhen  ^  is  approximately 
equal  to 


{-.^M 


lo  zmtn^i  r.Z 


•10 


(a. 5) 
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h     - 


it  n 


where  n  and  m  ^^  0  are  integers  such  that  0  <  -w—  <  1. 

The  maximum  departure  of  the  stream  lines  from  y  =  hA.  occurs 
at  X  =  0,   When  a  is  between  0  and  it  the  maxiinum  departure  occurs  at 
the  free  surface  and  the  behavior  is  like  the  ordinary  solitary  wave 
in  a  medium  of  constant  density.   Otherwise,  the  extrema  of  the  stream 
lines  with  respect  to  y  =  h'L.  occur  at  x  =  0,  and  the  values  of 
which  satisfy 


k  sin  a'L+a  cos  a/4,  =  0  , 


or  (since  ■-  =  tan  ^) 


cos  a(^L-i)  =  0 


W 


e  denote  these  values  by  't . .   They  are  approximately  equal  to 

A      -  2i+l 
•■1    ijja 

where  m  is  a  positive  integer  and  i  =  0,1,2, ... ,2m-l.   The  maxima 
and  minima  of  the  stream  lines  with  respect  to  y  =  h^t  are  given 
approximately  by 

^\  ^^^^    i 

qhe    sin  a^  =  qhe       (-1)   . 

The  extremura  nearest  the  bottom  is  a  positive  maximum;  and  the 
extrema,  v;ith  slightly  increasing  absolute  values,  alternate  between 
positive  and  negative  values  as  we  ascend  from  the  bottom  of  the 
channel  to  the  top.   For  obvious  physical  reasons,  an  extremum  must 
be  less  in  absolute  value  than  the  distance  from  its  place  of  occur- 
rence, (0,^.),  to  the  nearest  stream,  line  parallel  to  the  bottom. 


,>o   'to  .-r: 
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Hence  I'jhen  a  is  lar^'e  the  araplitudes  of  the  Internal  waves  and  the 
surface  wave  must  be  small.   Approximately  speaking,  the  stream  lines 
in  the  strip  0  <  ^  <  ^r-  are  lines  of  elevation  above  y  =  h^u   The 
lines  in  the  next  strip,  T=r-  <  'L  <  2*-^   are  lines  of  depression  below 
y  =  h^,    and  so  on  to  the  free  surface.   The  profile  of  the  free  sur- 
face is  given  by 


y  =  h+qhe   sin  a  sech'  "^wli-j^^ 

The  free  surface  is  a  wave  of  elevation,  but  if  a  is  not  between  0 
and  %   the  amplitude  tends  to  be  quite  small  since  sin  a  =  sin  2«mit  =  0, 
Let  us  denote  the  least  positive  root  of 

tan  a  =  — 
a 

by  a^.  For  k  small,  a,  is  given  approximately  by 

„2  AJ  k^ 

ex  -  2k  -  -V  . 

Per  this  root  the  stream  lines  are  lines  of  elevation  above  y  =  h^, 
and  the  maximum  araplitude  occurs  at  the  free  surface.   The  maximum 
amplitude  is 

q-n  =  qhe   sin  cu,  =  qhe   sin  ./Zk     . 

A  numerical  calculation  in  '.jhich  we  ultimately  retain  only  the  first 
power  of  k,  shows  from  (5.6)  that  the  wave  speed  is  approximately 
given  by 

(5.7)  |^^_l^(l-.$:)  . 

gh    ^^ 


A  similar  calculation  for  -m-|/3m  shows 


e&n 


-  j.>  /. 


J.?   !:'-cf 


odf!   i: 


■  w 
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Hence  the  stream  lines  for  a  =  a,  are  defined  by 


/t  A^     T'        ^  -1  v,^4.       k(>b-l)    sin  y2k^6  .    ^^  ,  2    x     /    ^  /^^  °- 
(5.8)      f(.x,Y)    =  h1+q,  e  ^^^ •    sech     -^  ^ 


sin/sic  ^^       he^  sin  /2k 

Notice   that    if  vie  take   the   limit   of  the  right  hand  side    of   (5.8)    as 
k  — >  0  we  have 

r(x,Y)   =  h'L+q]_>l  sech-  -^J-rr 

which  coincides  with  the  known  result  for  the  stream  lines  in  a 

medium  of  constant  density.   Also,  as  k  -*'  0,  the  formula  (5.7)  becomes 

c^  ^  ^1 

g-i       h 

which  is  the  Scott  Russell  formula. 

Case  II.   a  =  (2m+l)rL;  but  a  >  (2:i+1)7t:  (m  =  0,  1,  2...)  . 
For  a  value  of  a  which  is  nearly  an  odd  multiple  of  it  we  have 

cot^^-l   j,2.^2 
cos  a  =      -    =  -.^-^ 
cot  -^  + 1   k  +a 


and 


^  =     (a^+k^-<-2k)(a^-k^)(k^+9a^) 

^1    2a^  (  a'+k'^T78'( k^+a^ ) e'^  cos  a-(k^-a2)] 

=  -  (a^+k^+2k)(k--!-9a^) 
2a^(a^+k^)(8?^^oT' 


which  is  negative.   If  we  vjrite 


rri  ="  (y- 


i   3-.;  0  '.J-  ^ 


ssbiDntoo 


rjr;B 
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3m, 


the  stream  lines  are  given  by 


(5.9)     y  =  r(x,Y)  =  hA,-QJie^''^  sin  aA  sech^  ik^'^^^^^ 


where 


"^1  ^ I;ka^(k^+a^)-[8e^-l] 


In  terms  of  Q,  a,  and  k  the  speed  of  the  wave  is  given  by 

rr 

35 


^  =  ■f  -  5^  '^ 


or 


(5.10)  S|  =  £!^  Cl  .   |I:a3f8e^-l)q 

c^     "-'^       (a''+k'^+2k)(k'^+9a'^) 

From  (5.9)  we  see  that  the  stream  lines  are  parallel  to  the  horizontal 
bottom  of  the  channel  for  those  values  of  -^  x;hich  satisfy  sin  a^-  =  0. 
These  horizontal  stream  lines  occur  when  ^ is  approximately  equal  to 


2in+l 


where  n  takes  the  values  n  =  0,1,2, ••• ,2m+l. 

The  extrema  of  the  stream  lines  with  respect  to  y  =  hK.  occur  at 
X  =  0  and  those  values  of  /t.  vjhich  satisfy 

k  sin  o/T+a  cos  a^  =  0  . 

Since  for  the  case  under  consideration  -cot  %  ~  ■„   t    these  values  of 
*l.    satisfy 


Sii 


2    Oi:.' 


. /c^oo 


«  rt 


ds   'iiLfoco  .,  rl  =   v. 


7k 


sin  a{K-  i)  =  0  . 


Vie  will  denote  these  values  by'L.   They  are  approximately  equal  to 

-  2J+1 


K  zz 


2T2rQ+l ) 


when  j=0,l,2, . , ,  ,2ra.   The  maxima  and  minima  of  the  stream  lines  with 
respect  to  j-  =  h<L  are  given  approximately  by 

-Qhe^  J  sina^u  =  -(-l)^  Qlie^^SS^^  ^ 

J 

The  extreraum  nearest  the  bottom  is  a  negative  minimum;  and  the 
extreraa  with  slightly  increasing  absolute  values  alternate  between 
negative  and  positive  values  as  we  ascend  from  the  bottom  of  the 
channel  to  the  top.   The  absolute  value  of  an  extremum  must  be  less 
than  the  distance  from  (O,'^;,)  to  the  nearest  horizontal  stream  line. 
This  physical  condition  imposes  an  upper  bound  for  Q.   The  amplitude 
of  the  internal  waves  must  tend  to  zero  as  a  increases. 

For  the  strip  approximately  given  by  0  <  ^  <  p— rr  >  the  stream 
lines  are  lines  of  depression  below  y  =  h'L*   In  the  next  strip 

2m+T  "^  ^  "^   "2m+l  ^^®  stream  lines  are  lines  of  elevation  above 
y  =  h^  The  alternation  of  strips  of  lines  of  depression  and  strips 
of  lines  of  elevation  continues  until  the  free  surface  is  included. 
The  profile  of  the  wave  at  the  free  surface  is  given  by 


y  =  h-Qhe^  sin  a  sech^  ^E^'sm^^   • 

This  is  a  wave  of  elevation  because  a  is  slightly  larger  here  than  an 
odd  multiple  of  %   and  hence  sin  a  is  negative,  but  since  sin  a  is 
small  the  amplitude  of  the  wave  tends  to  be  small. 


acv; 


'.     fid     CJtiiWili    i* 


,  2Xfo     ' 


— eri    '5* 
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As  a  specific  example  of  a  solitary  wave  of  the  kind  v/e  are  now 
considering,  let  us  take  the  solitary  wave  which  corresponds  to  the 
least  value  of  a  which  satisfies 

-cot  2  -  a  • 
We  will  denote  this  value  by  a^,   VJhen  k  is  small,  ttp  is  approximately 

r,     <^'  ^  J.  2k 
^2  -  ^  "■  If  • 

m, 
For  this  value  of  a,  a  calculation  of  -7^7=-,  in  i.'hich  we  retain  only 

the  first  power  of  k,  gives 


^^1  _  28kTi 


o 

and  the  stream  lines  are  given  by 


y  =  r(x,Y)  =  h  -Qhe^'^-  sin  a.^t  sech^  ^/ogSkir   ^ 
We  see  that  there  is  a  horizontal  stream  line  when  ^  is  such  that 


sin  <iJl   =  0  . 


This  value  of  ^l   is  approximately 

/-.  _     It    '^      %        2J  T  2k 

^      a^       rr2K       -^  "  ~2     • 

2         7t+— —  Tt 

II 

Thus  the  horizontal  stream  line  is  close  to  the  free  surface.   In 

the  strip  0  <  ''l  <  77-  the  stream  lines  are  lines  of  depression  below 

y  =  hA.   In  the  narrow  strip  — -  <  'Z.  <  1,  the  stream  lines  are  lines 

^2 
of  elevation  above  y  =  h^.   The  stream  line  of  greatest  depression 

below  y  =  h^J.  is  given  by  the  value  of  /l  which  satisfies 


■+ ,- 
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sin  a^(^L-i)  =  0  , 


namely  4  =  l/2.  The  rnaximum  departure  of  this  streara  line  from  y  =  -^ 
is  Qhe  '  sin  -w-  =  0-  say.  The  maximum  amplitude  of  the  vrave  at  the 
free  surface  is 

k  2     °-Z 

-Qhe  sin  a^  =  -Q-  2e  cos  -^   . 

k 

This  is  approximately  2Qp  e   sin  —  which  is  small  compared  with  Qp. 

The  approximate  speed  of  the  solitary  wave  can  be  obtained  from 

(5.10).   In  terms  of  Qp  x;e  have 

c      "^^^        (a^+k''+2k)(k'^+9a^)     §     a^ 

"^  ^    he"^  sin  -^ 

If  we  perform  the  calculation  and  retain  only  the  first  powers  of 

k  we  find 

c^  'jy  2k  r,  .  28k  ^2i 
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